Trace formulae for non-equilibrium Casimir interactions, heat radiation and heat 

transfer for arbitrary objects 



Matthias Kriiger, 1 Giuseppe Bimonte, 2 Thorsten Emig, 3 and Mehran Kardar 1 

1 Massachusetts Institute of Technology, Department of Physics, Cambridge, Massachusetts 02139, USA 
2 Dipartimento di Scienze Fisiche, Universita di Napoli Federico II, 
Complesso Universitario MSA, Via Cintia, 1-80126 Napoli, 

Italy and INFN Sezione di Napoli, 1-80126 Napoli, Italy 
3 Laboratoire de Physique Theorique et Modeles Statistiques, 
CNRS UMR 8626, Bat. 100, Universite Paris-Sud, 91405 Orsay cedex, France 

(Dated: July 17, 2012) 

We present a detailed derivation of heat radiation, heat transfer and (Casimir) interactions for N 
arbitrary objects in the framework of fiuctuational electrodynamics in thermal non-equilibrium. The 
results can be expressed as basis-independent trace formulae in terms of the scattering operators of 
the individual objects. We prove that heat radiation of a single object is positive, and that heat 
transfer (for two arbitrary passive objects) is from the hotter to a colder body. The heat transferred 
is also symmetric, exactly reversed if the two temperatures are exchanged. Introducing partial 
wave-expansions, we transform the results for radiation, transfer and forces into traces of matrices 
that can be evaluated in any basis, analogous to the equilibrium Casimir force. The method is 
illustrated by (re)deriving the heat radiation of a plate, a sphere and a cylinder. We analyze the 
radiation of a sphere for different materials, emphasizing that a simplification often employed for 
metallic nano-spheres is typically invalid. We derive asymptotic formulae for heat transfer and non- 
equilibrium interactions for the cases of a sphere in front a plate and for two spheres, extending 
previous results. As an example, we show that a hot nano-sphere can levitate above a plate with the 
repulsive non-equilibrium force overcoming gravity - an effect that is not due to radiation pressure. 

PACS numbers: 12.20.-m, 44.40. +a, 05.70.Ln 



I. INTRODUCTION 

Quantum thermal fluctuations of electromagnetic 
waves lie at the heart of statistical physics, accounting 
for seminal discoveries such as Planck's law for thermal 
radiation of black bodies, introduced almost a century 
a §° [H- The equilibrium Casimir force 0] between par- 
allel metallic plates can be equivalcntly attributed to the 
fluctuations of the electromagnetic field, or to the charge 
and current fluctuations in the plates @. At separa- 
tions much smaller than the thermal wavelength, which 
is roughly 8 yum at room temperature, these forces are 
generally dominated by quantum zero point fluctuations, 
whereas at larger separations thermal fluctuations also 
need to be considered 0, 0] ■ I n thermal equilibrium the 
tools of statistical physics can be exploited to ascribe a 
(Hclmholtz) free energy to a collection of objects; deriva- 
tive of the free energy with respect to separation (or ori- 
entation) of the objects yields the force (or torque). The 
electromagnetic free energy can itself be compactly ex- 
pressed in terms of the scattering operators of the objects 
(see, e.g., Refs. [l-Q). 

A notable property of the equilibrium formulae is that 
integrals over frequency can be evaluated along the imagi- 
nary axis [|| , where expressions for the response functions 
are much smoother, and forces are not very much influ- 
enced by the precise position of material resonances (in 
contrast to thermal non-equilibrium) . Another feature of 
equilibrium forces is that stable situations are impossible 
under rather general conditions, e.g., the free energy as 



a function of the assembly of objects in vacuum has no 
minima Q. (Unstable repulsion can still be obtained in 
certain cases 0].) 

The improved precision of measurements of force and 
heat transfer at sub-micron scales have provided renewed 
incentive in the past decade to examine fiuctuational 
electrodynamics (FE) for objects at different tempera- 
tures [l(J, [H| (but also other forms of non-equilibrium, 
e.g., objects in motion [12]). The current status of the- 
ories for thermal non-equilibrium have in general two 
commonalities: additional assumptions about the system 
have to be made (e.g. the assumption of local equilibrium 
within each object [l3j ) to enable any prediction, and 
the basic tools of statistical physics, e.g., the free energy, 
cannot be employed complicating analysis. Assuming lo- 
cal equilibria, the current fluctuations in each body are 
treated separately at the objects' temperature, e.g., using 
fiuctuational electrodynamics introduced by Rytov over 
60 years ago [l^. Recently, out of equilibrium Casimir 
forces have been computed in a number of cases includ- 



ing parallel plates [ll(, deformed plates 14 L as well as 



a plate and an atom in different setups 15-tl7l|. There 
also exists a large body of work on forces between atoms 
or molecules in non-equilibrium fl8l - |2"4"j ] . Formalisms for 
treating FE for arbitrary objects at different tempera- 
tures have been recently presented [25t - [28l | . In partic- 
ular, for compact objects, radiation from the environ- 
ment contributes to the force and has to be incorporated. 
Quite generally thermal non-equilibrium can be repul- 
sive [H, H [H, HlHil. and allow for stable zero force 
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points [27|, l29( . Repulsion can even occur at separations 
far below the thermal wavelength if the resonances of the 
materials are suitably detuned [111 . It is also possible for 
a hot and cold sphere to exert equal mutual forces in the 
same direction, leading to a self-propelled state [27| • Two 
parallel nanotubes where shown to be good candidates 
for experimental detection of non-equilibrium forces, as 
these can be made relatively strong [30(. Recently, the 
high temperature limit was investigated in thermal non- 
equilibrium (3l| . 

Heat radiation and transfer are of particular interest 
when the size and/or separation of the objects is com- 
parable to, or smaller than, the thermal wavelength, be- 
cause then they differ strongly from the predictions of the 
Stefan-Boltzmann law. For example, as confirmed exper- 
imentally f32l - l3~H , there is a considerably larger near-field 
heat transfer due to tunneling of evanescent waves. The- 
oretical computations of heat transfer were only recently 
extended from two parallel plates [HJ or dipoles [36[ 
to compact objects of finite size, more precisely to two 
spheres [37], HH and a sphere in front a plate 2g, [H, H(| ■ 
Numerical studies for objects whose scattering properties 
are not known analytically include transfer for periodic 
structures [4l[ , as well as a cone or finite cylinder in front 
a plate , and very recently, numerical scattering tech- 
niques were implemented more generally (42j . (See also 
Refs. [43l - [49| for recent studies of various aspects of heat 
transfer.) 

The radiation of single spheres and plates has theo- 
retically been studied by many authors [2(| I50l - l55| . 
Radiation of single cylinders, for which an early calcula- 
tion by Rytov exists [5(j, was only recently formulated 
in terms of scattering theor y 12(1 Il54l [H[ , and has been 
investigated experimentally [56|-|60(, mostly focussing on 
polarization effects in the emission of thin tubes (cylin- 
ders). Recently, rotating objects were considered, which 
emit spontaneously at zero temperature [6l|, [62j . Refer- 
ence [63| provides a basis independent trace formula for 
the radiation of an isolated object (see also [62|). 

The main result of this paper is the derivation of gen- 
eral trace formulae for heat radiation, transfer and non- 
equilibrium forces (including contributions from the en- 
vironment) for arbitrary objects. The trace formulae do 
not refer to any particular (wave) basis and hence can be 
employed in rather general situations. We give proofs for 
the positivity of radiation and transfer. To demonstrate 
the power of our general results, we provide analytical as 
well as numerical examples. The paper is divided into 
three main parts: The first part (Sees. IIHI - IVI|) formu- 
lates the problem in terms of scattering operators (e.g. 
T(r, r')). In the second part (Sec. IVII|) . we transform 
these expressions into matrix forms in arbitrary partial 
wave basis. The third part (Sees. IPXl - IXII[) is devoted to 
specific analytical and numerical examples. 

In particular, in Sec. [II] we describe the model and de- 
rive the non-equilibrium correlation function of the elec- 
tric field, followed by an introduction of the T opera- 
tor, in sec. nni we derive formulae for both the elec- 



tric field correlation as well as the emitted energy for a 
single object in terms of T, and prove that the emitted 
energy is positive for any object made of passive ma- 
terial. In Sec. IIV1 we provide the non-equilibrium field 
correlator for N objects. This correlator is then used to 
derive a trace formula for heat transfer in Sec. |Vj where 
we also prove the positivity of the transfer as well as 
its symmetry with respect to a permutation of tempera- 
tures. Trace formulae for the Casimir force are derived 
in Sec. EH In Sec. I VIII we introduce partial wave ex- 
pansions and express the trace formulae in terms of the 
corresponding matrix expressions. A short discussion of 
differences between equilibrium and non-equilibrium cal- 
culations follows in Sec. IVIIII In Sec. [DO we (rc)derive 
the radiation of a plate, a sphere and a cylinder and an- 
alyze the radiation of a sphere in detail. Then we give 
asymptotic expansions for the cases of two spheres and a 
sphere in front of a plate for heat transfer (Sec. |X| an d 
forces (Sec. IXlj) . Sec. IXIII provides numerical examples 
for the sphere-plate case, demonstrating that the non- 
equilibrium force can lead to stable levitation points. In 
App. [A] we show that the results for equilibrium forces 
can be derived from our trace formula for non-equilibrium 
forces. Appendices [BlIEl present the partial wave expan- 
sions and conversion matrices used in Sees. IIXHXII 

We note that during submission of this manuscript, 
we became aware of an independent, but partly related 
work on a trace formula for heat transfer, its symmetry 
and positivity (42j. 



II. NON-EQUILIBRIUM FLUCTUATIONS 
GENERAL CONCEPTS 



We consider an arrangement of N objects labelled by 
a = 1 . . . N , in vacuum at time- independent, homoge- 
neous temperatures {T a }, and embedded in an environ- 
ment at temperature T cnv . The objects are characterized 
by their electric and magnetic response e(w;r, r') and 
p(lo; r, r'), which can in general be nonlocal complex ten- 
sors, c(uj; r, r') = Eij(u); r, r'), depending on frequency uj. 
In this non-equilibrium stationary state, each object is 
assumed to be at local equilibrium, such that the cur- 
rent fluctuations within the object satisfy the fluctuation 
dissipation theorem at the appropriate temperature [l3j . 

For any two field operators A and B, we consider the 
symmetrized expectation value, 



(A(t, r) B(t',r')) s = l -(A{t, r)B(t\ r') + B(t\ v')A(t, r 

Symmetrization ensures the reality of (A(t, r)B(t', r')) s 
for generally non-commuting quantum operators. In sta- 
tionary conditions, the expectation value depends only 
on the time difference t — t', and we can define the spec- 
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tral density (A(r)B* (r')) u by 



A(t,T)B(t',r')) = 



dio 
2tt 



-Mt-t') ^4( r )£*( r ')) D 



(2) 

Due to the reality of the correlation function on the left 
hand side, the real part of the spectral density is an even 
function of frequency, while its imaginary part is odd. 

The relevant quantity for our considerations is the 
spectral density C = Cy of the electric field E at points r 
and r', from which heat radiation and Casimir forces can 
then be extracted (see Sees. IVl and I VII below) . defined by 



C<,(r,rOM^(r)£*(r')>o 



(3) 



In the following we shall not make explicit the depen- 
dence on u) in frequency dependent quantities. In order to 
derive this correlation in the considered non-equilibrium 
state, we start with the equilibrium case where all tem- 
peratures are equal, T a = T env = T. Then, C is well 
known, and can be expressed in terms of Gjj , the dyadic 
retarded Green's function of the system 64|. This rela- 
tion is a variant of the fluctuation dissipation theorem, 
and reads 1811 



C%(T; r, r') = [a(T) + a ] Im G« (r, r'). 



(4) 



Here 



a(T) ee sg nH^![exp(%|/fc s T) - l]" 1 , (5) 
cr 

contains the occupation number of modes with frequency 
ui, c is the speed of light, and H is Planck's constant. Zero 

point fluctuations have amplitude ao = sgn(a;) 4 ^" , but 
play no role in non-equilibrium phenomena- they are 
independent of temperature and present everywhere in 
space. They do contribute to the total force, as they are 
also responsible for the zero point Casimir effect. 

The dyadic Green's function obeys the Hclmholtz 
equation 



fr,r') = M(%r- 



where 



= V x Vx, and \8c 



V=^(c-I) + Vx [I-- 
c 2 \ P 



(r-r'), 



Vx 



(6) 



(7) 



is the potential introduced by the objects. We can sepa- 
rate the term (E ® E*)° = ao ImG (where we introduce 
the dyadic vector notation (E <g> E*)^- = E t E*) involv- 
ing the zero point contribution, and concentrate on the 
remaining T-dependent terms. These can be _split up ac- 
cording to their originating thermal sources [26[ , yielding 
N + I terms, including the contribution from the envi- 
ronment. The equilibrium correlation in Eq. Q can thus 
be rewritten as 

C eq = (E g> E*)° + J2 C a( T ) + C onv (T). (8) 



This form shows the different contributions to the electric 
field correlations in equilibrium. The term C* C (T) repre- 
sents the radiation from the sources in object a, and is 
given by 



C* C (T) = a(T)GIm[V Q ]G*. 



(9) 



V a is the potential of object a, i.e., e and p in Eq. (J7J arc 
replaced by <e q and p a . Equation (j9|) is identified with the 
original definition of the field correlator in Rytov's for- 
malism [H, H3l- For the case of Im[[Ua] = 0, we have the 
more familiar expression C s a c (T) = o(T)^GIm[e a ]G*. 
One difference to the heat radiation of object a as de- 
scribed in Rcfs. [l3|, EH is that in Eq. ([9]) the radiation 
is scattered by all objects, such that Im[V a ] is multiplied 
from both sides by the full Green's function G. For a sin- 
gle object in isolation [see Eq. (|22]l below], G is replaced 
by G a , the Green's function of object a in isolation. 

Comparing Eqs. ^ and ((5J), we derive the last term 
in Eq. ([8]) which can be identified as the contribution of 
the sources in the environment, 



C cnv (T) = -a(T)GIm [G^ 1 ] 
This can be evaluated further by using 



Im[V] 



Imf 1 



(10) 



(11) 



which follows directly from Eq. where Go is the free 
Green's function solving the wave equation for V = 0. 
Fluctuations in the vacuum are taken into account by 
the nontrivial term Gq , which can be attributed to in- 
finitesimal environmental "dust" (64J. Integrating such 
infinitesimal "dust" sources over the infinite space of the 
environment yields a finite result. 

The intuitive result that the equilibrium field in Eq. © 
is the sum of the radiation emitted by the sources in the 
objects and the sources in the environment can also be 
corroborated by deriving C env in Eq. (JSj) along another, 
straight forward route. We start from the field E sourced 
by the environment, without any objects present. It has 
the correlator 



E®E*)£ ee = a(T) ImGc 



(12) 



Then, adding the cold, i.e., non-radiating objects, 
generates scattered fields according to the Lippmann- 
Schwinger equation. If the field E solves the Hclmholtz 
equation in free space, then the following E sc solves it 
with the objects present, 



(13) 



Applying the operator GG 1 to both fields in Eq. (I12p , 
one finds for the environment contribution 



C cnv (T) = 

= -a(T)GIm [Gq 1 ] G*, 

reproducing the last term of Eq. (|5J) . 



(14) 
(15) 
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Having identified the contributions from the different 
sources, we can now change the temperature of these 
sources independently, denoting by T a the temperature 
of object a and by T env the temperature of the environ- 
ment. The field correlator in the non-equilibrium situa- 
tion is then a simple modification of Eq. (JSj> to 



C ne< i({T a },T env ) = (E®E*)° +J2<C s a c (T a ) 
+ C cnv (T cnv ). 



(16) 



Equation (fl6|) gives the general field correlator for arbi- 
trary combinations of the temperatures of the N body 
system. It contains N + 1 unknown terms, due to the 
N + 1 sources in the system. We note however that one 
of the sources can be eliminated by introducing the equi- 
librium correlation function at finite temperature. We 
chose to eliminate the environment contribution and ob- 
tain m 

C ncq ({T Q },T cnv ) = C cq (T cnv )+]T [C s a c (T a ) - C s a c (T cnv )} . 

a 

(17) 

This form shows that we only have to evaluate the N 
terms C* c (assuming C oq is known), in order to compute 
heat transfer and forces depending on TV + 1 tempera- 
tures. For Casimir forces, C eq (T env ) will give the equi- 
librium force at temperature T cnv . However, due to its 
equilibrium origin, this term does not contribute to the 
heat transfer. In Sec. IIV1 we give the final formula for 
C neq ({T a }, T env ) in terms of the scattering operators of 
the objects. In Sees. [V] and IVIi we shall derive the re- 
sulting heat transfer and forces, respectively. 

We conclude this section by introducing the classical 
T-operator which provides a convenient way of rewriting 
the Hclmholtz equation as a Lippmann-Schwinger equa- 
tion (65| . Starting from 



E sc = E + G VE S 



(18) 



we can formally write E sc in terms of the T-operator [84 
as 



E sc = E + G TE. 
Solving for T, we obtain 

1 

T = V- 



(19) 



(20) 



which to the lowest order equals V, as in the Born ap- 
proximation. Comparing Eq. (|T9l) to Eq. (fl"3)) . we find 
the following relation between T and G Q, 



G = G + G TG . 



(21) 



Note that T is the scattering operator of the entire col- 
lection of objects, whereas we shall use T a for object a 
in isolation. 



III. RADIATION OF ONE OBJECT IN 
ISOLATION 

In this section we derive the correlator C Q of a single 
warm object in a cold environment, a prerequisite for 
the derivation of C^ c needed in Eq. (|17|). We shall also 
compute the energy emitted of the isolated object and 
prove that it is positive. 



A. Field correlations 

Equation (fT7| requires the expression for C^ c , the field 
sourced by object a and scattered by all objects. In the 
absence of other objects, the thermal field correlator (ne- 
glecting zero point terms) satisfies 



C a (T a ) = a(T a )G a lm[V 



(22) 



C Q (T Q ) differing from C^ c in Eq. © by the appearance 
of the Green's function G Q instead of G. There are dif- 
ferent ways to evaluate Eq. (f2"2"|) for a specific geometry. 
The most straightforward approach is to start with the 
Green's function G a with one point inside the object and 
one point outside, as this is the structure of Eq. (|22[) : 
Im[V Q ] is only nonzero inside the object and we are in- 
terested in the field outside. Once this is accomplished, 
Eq. (|2"2"j) can be directly evaluated by an integration over 
the volume of the object. Such an approach was used, 
e.g., in Refs. [Til |37| to find the non-equilibrium Casimir 
force for parallel plates and the heat transfer between 
two spheres, respectively. 

Since we aim to describe the non-equilibrium effects 
through the scattering formalism, we would like to ex- 
press desired observables - starting from the heat radia- 
tion of a single object - in terms of the scattering oper- 
ators {T Q }. The single object's radiation can indeed be 
expressed in terms of T a by starting from the equilibrium 
situation of T a = T onv , where the field correlator C^(T a ) 
(not containing zero point fluctuations) can be split [in 
a manner similar to Eq. ([5])] into contributions of heat 
sources from the object and from the environment, as 

= C a (T a ) + C™ v (T a ) = a(T a ) Im G Q , (23) 
C« nv (T Q ) = a(T a )G a Im [Gq 1 ] G* . (24) 

The last equality in Eq. (l23l) follows from Eq. ((4]) after 
reduction to one object [64| . In order to arrive at the 
desired heat radiation of the object, we solve Eq. (|23| for 
C a (T a ), yielding 



c a (r a ) = -cr(T a ) + TO, 



(25) 



Cq 9 (Tq,) is readily expressed in terms of the T-operator 
via Eq. (f2"Tj) . It remains to express the radiation sourced 
by the environment, see Eq. (|24|) . in terms of T. This 
can be achieved along different routes. We want to first 
present the way introduced in Ref. [26| and then present 
a general formula for the radiated field. 
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1. Integration over environment "dust" 

Following the interpretation introduced in Ref. [64 |, 
the environment can be regarded as composed of "dust" 
characterized by a homogeneous dielectric response e env 
in the infinite space complimentary to the object. The 
Green's function of the system [for V = ^-I(e en v — 1) 
outside, and V = V a inside the object] is denoted G Q . It 
is a simple modification of G Q as a uniform e cnv — 1 only 
changes the speed of light outside the object, replacing c 
with c/ ^/e cnv . Formally interpreting the environment as 
an additional object, we can use Eq. © to get the fields 
sourced by it. Taking e onv — > 1 yields a well defined 
radiation from the environment, 



-T(T a ) = lim , a(T a )—G a Im[ £env ]G* 



(26) 



Writing out this equation explicitly, it is possible to see 
more clearly the operator products involved, as 

C°%-(T Q ;r,r')-a(r Q )^ 

x lim V/ d 3 r"G Q , jfc (r,r")Im[e cnv ]G; ifej (r",r'). 

(27) 

First, we note that all the arguments of the Green's func- 
tions in Eq. (|27|) lie outside the object, such that G Q can 
be found by use of scattering theory (see Sec. IVIIj) . One 
additional simplification occurs since the limit e cnv — > 1 
allows to neglect any finite region of integration. We can 
hence restrict the integration range to the region with 
£i(r") > max[£i(r),£i(r')], where £i is the radial com- 
ponent which distinguishes the two expansions of Go in 
Eq. ([79")) below. This practical simplification which holds 
for any finite r, r', allows restriction to one of the cases 
in Eq. (|T9"|) . such that finally 

C^(T a -v,r') = s lim i J2 J d 3 r" 

k f(r»)>max{£(r),£( r ')} 

G a , lfe (r, r") Im[e env }G*^ kj (r", r'). (28) 



Equation (f2"5)l was presented in Ref. [26( , and its applica- 
tion for a cylindrical object was worked out in detail in 
Ref. [Ill . It can be evaluated in a straightforward manner 
if the matrix elements 7^' [see Eq. (|82p] are known. 



2. A general formula for the field correlations 

A simpler method for expressin g t he radiation of 
the environment, presented in Ref. [55|, allows to give 
Eq. (|25[) in closed form. First, we rewrite the expression 
for the environment radiation, Eq. (fiM)) , using Eq. (f2~Tj). 
and (after a few steps) find [55[ 

C c a m (T a ) = a(T a )[l + G T Q ]ImGo[T*GS + 1]. (29) 



Here, we have used that the T operator is symmetric [be- 
cause of Eq. (f2"Tj) and since G [64[ is symmetric] . Using 
this form, and additionally writing Cq? in terms of Tq, 
via Eq. (f2~Tj) , the heat radiation in Eq. (|25p can be given 
in closed form for an arbitrary object. After some ma- 
nipulations wc find, 



t (T a ) = a(T a )G 



(T* -T a )-T a Im[G ]T* 



(30) 



This describes the radiation of an arbitrary object in 
a basis-independent representation in terms of two well 
known quantities, the free Green's function and the 
T-operator. Such a basis-independent representation 
can be of advantage when numerical methods are em- 
ployed [4(j to find T a . We emphasize again that it holds 
for any material properties, with magnetic or electric 
losses. In its derivation, which requires only Eqs. (|23[) 
and (|29|) . we even do not have to be cautious about (mag- 
netic or electric) material losses as this information is 
contained in the T-operator. In case of lossless materi- 
als, Eq. (ptUf does not correspond to a field configuration 
that supports energy transport, as will be demonstrated 
below. To simplify notation, we define the radiation op- 
erator R„ and write 



C Q (T Q ) = a(T a ) 



with 



i a — Go 



Im[T Q ] - T Q Im[G ]T* 



G,"i 



(31) 



B. Trace formula for the emitted energy 

In the previous subsection we derived the correlation 
function of the electric field for a single object in an envi- 
ronment at zero temperature, described by the radiation 
operator R Q in Eq. (f5U|) . We now consider the heat emit- 
ted by this object, obtained by integrating the normal 
component of the Poynting vector S over a surface E a 
enclosing the object, as 



(32) 



(33) 



H a = J S • n Q , 

with 



It is straightforward to evaluate the integral in Eq. ([32]) 
for simple objects such as a sphere or a plane. Never- 
theless, a more convenient and illuminating form can be 
achieved by reconsidering the derivation of the Poynting 
theorem [661 ] , starting from the work done by the fluctu- 
ating fields in a volume element located at r. This work is 
given by the electric field times the total electric current 
J at r. Hence the total work done on the object (where 
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we include a minus sign to get the emitted energy) is 
given by an integral over the volume of the object, 

/OO 1 P 
— / <?r(E(r).J*(r)) u . (34) 
- oo 2lT JV a 

To proceed, we use the free Green's function to convert 
between the total field and the total current, 

E = 47ri^-G J. (35) 

Applying this to Eq. ([34]) . leads to 

! £36) 



The heat transfer _£f from the object to the environment 
at finite temperature T onv follows directly from detailed 
balance, as 

H = H a (T a ) — H a (T onv ) . (38) 

Note that the trace in Eq. (f3"7| of the operator IRcGq -1 * = 
(R q Gq 1 *)i J -(r,r') is both over the vector indices i and 
j, as well as the positions r and r' [85| . This can be 
converted into a more familiar trace in a partial wave 
basis, see Eq. (p?5j) in Sec. IVIII below. However, when 
Im[Go] is expanded in wave functions it contains only 
propagating waves, such that the trace in Eq. (|95|) is 
restricted accordingly. 

C. Positivity of heat radiation 

Causality implies that the potential V(w) in Eq. Q 
is an analytic function in the upper half of the complex 
frequency plane, with the property 

W(-z*)=Y*{z). (39) 

For real frequencies, the above condition implies that the 
real part of V(w) is an even function of uj, while its imag- 
inary part is odd. By virtue of this symmetry, we may 
consider only positive frequencies. The imaginary part of 
the potential V(w) for a body made of a passive material 
must be positive semi-definite, i.e. 

Im[V] > . (40) 

For any positive semi-definite operator A, the product 
BAB^ is also positive semi-definite, and we have 

Go *G Im[V] G*G^ 1 > . (41) 



From Eq. (|30|) we note that the correlation function of the 
electric field carries the Green's function Gq on its right- 
most position. Hence the operation of Gq 1 * in Eq. (|36|) is 
easily performed, leaving T a on the rightmost position. 
We now note that the range of integration can be ex- 
tended to all space, as MaGg -1 * is nonzero within the vol- 
ume V a only. This is because T Q (r,r') is only nonzero if 
both arguments are located within the volume V a . Thus, 
the integral in Eq. (|36|) together with the scalar product 
turns into a trace of the operator M a Gg 1 *, and we get for 
the energy emitted by an arbitrary object with scattering 
operator T a , 



(37) 

I 

If we furthermore use Eq. ([20]) , we directly find 

Im[T] -TIm[G ]T* > . (42) 

As Im[Go] is a positive semi-definite hcrmitian operator, 
and the product of two positive semi-definite hcrmitian 
operators is also positive semi-definite, the operator to 
be traced in Eq. (|37l) is positive semi-definite as well. 
This shows that the emitted energy H a is a nonnegative 
number, 

H a >0. (43) 

While this is expected on physical grounds, to our best 
knowledge, it has not been proven for an object of arbi- 
trary shape before. Furthermore, as Tlm[Go]T* > 0, we 
conclude from Eq. (|4*2"j) 

Im[T] > TIm[G ]T* > , (44) 

which proves that Im[T] is a positive semi-definite oper- 
ator. 



IV. RADIATION FROM MULTIPLE OBJECTS 
AT DIFFERENT TEMPERATURES 

In Sec. piip we derived C Q (T a ) for the field radiated by 
an isolated object a, in terms of the radiation operator 
M Q . This radiation is scattered at all other objects in the 
system, leading to the modified correlator C Q C , which is 
the unknown term in the total non-equilibrium correlator 
of Eq. © . In the following, V a will denote the composite 
potential of all objects except object a. For two objects, 
V g is the potential of the second object. 

We place the cold objects described by V a into the field 
radiated by object a. If the solution to the Hclmholtz 



H a {T a ) = — / du^^ Tr {Im[G ] Im[T a ] - Im[G ]T a Im[G ]T a } . 
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equation for object a alone is denoted by ~E a .iso, then 
the solution E Q for all objects can be expressed through 
the Lippmann-Schwingcr equation as [65| 



E„ 



(45) 



Writing the potential V g in terms of the Green's function 
G a , we arrive at [26l 



E„ 



^oi^a.iso i 



(i + c T a ; 



i 



1 - G T Q G T e 



(46) 
(47) 



The multiple scattering operator O a depends on the com- 
posite T-operator T a describing scattering by the other 
objects, as well as on T Q . Expanding the denominator of 
0„ leads to 



(l + G T a )[l + G T Q G T e 



(48) 



where the terms in square brackets correspond to an in- 
creasing number of back and forth scatterings between 
the objects. The expansion of Eq. (|48|l can be useful in 
order to get simplified analytical results in certain cases, 
as will be shown below. We note, however, that this 
expansion does not necessarily converge at close sepa- 
rations, [8r| and a multiple scattering expansion might 
not be as useful as in equilibrium situations where a fast 
convergence of such series is observed |67j . 

Applying the multiple scattering operator to both sides 
of the field correlator describing the radiation emitted by 
object a, we arrive at the final formula for the correlator 
as 



C„ C (T Q ) = O a C a (T a )Qt = a (T a ) 



Q Q . (49) 



In the last line, we have used the definition of the radia- 
tion operator in Eq. pip . Equation (|49p. together with 
Eqs. ® and (T7|) yields 

C nct H{T Q },T onv ) = C«HT cnv )+^ [Cl c (T a ) - C s a c (T cm )} . 



This constitutes our final result for the field correlator 
in the considered non-equilibrium situation, expressed 
in a basis-independent representation in terms of the T- 
operators of the objects. 



V. TRACE FORMULA FOR HEAT TRANSFER 

A. Two objects 

Let us consider two objects labelled 1 and 2 at tem- 
peratures T\ and T2, respectively. One can define differ- 
ent energy fluxes, which in general also depend on the 
temperature of the environment. For example, an exper- 
imental setup could measure the total energy absorbed 
by object 2 close to object 1 in an environment at a yet 
different temperature T cnv . While we derive the result 



for all possible cases in Eq. (|67[) below, we first focus 
on the usual definition in literature [U Hil H3] , where 
T cnv is assumed to be zero (or irrelevant in the near field 
regime), and one considers only the energy exchanged be- 
tween two objects. The component of radiation emitted 
by object 1 and absorbed by object 2, will be indicated 

(2*) 

by heat transfer rate H\ (?i). In turn, the emission by 2 

which is partly absorbed by 1 is quantified by H^\t 2 ). 

In the previous section we derived the correlation func- 
tion of the electric field. For computing the transfer 

(2) 

rate H{ , the standard method is to integrate the nor- 
mal component of the Poynting vector [related to Cf c in 
Eq. (|49p ] over a surface E 2 enclosing only object 2. As 
in Sec. IIIIBI for the heat emitted by a single object, we 
prefer to recast the transfer in terms of a volume inte- 
gral, which can then be turned into a trace. We have 
to evaluate again the expression (E(r) ■ [Gq *E*)(r)) . 
But now, aiming at the absorption by object 2, the in- 
tegral must be restricted to the volume of object 2, and 



the contribution from 



becomes 



7T ./ ok B T x 







dPr( 



^Go 1 *)^!-,!-) 



(50) 



We included a minus with respect to Eq. to indicate 
the energy absorbed by object 2. In the second line, we 
have replaced the correlator by Cf c , resulting in the op- 
erator OiRiO|Gq : *, which has to be traced over vector 
components i [because of the dot product in Eq. ([34]) ] 
and integrated over the volume V2. It can be split up as 



DiRiC^ 1 " 



with the two parts 



1 



1 - T*G^T*G5 

I , 

1 - T*G5T*G5 



T 



2 > 



(51) 

(52a) 
(52b) 



where i and s stand for "interaction" and "self," respec- 
tively. The splitting into the two operators in Eq. ([521 is 
done because they differ precisely by the operator on the 
most right position, 



= ...T 



2 > 



h^i®\ s = ...T* (or . . .Ti) , 



(53) 
(54) 



where the dots stand for the remaining parts of these 
terms. Now we note that O1M1OJ ^r, r) is only nonzero 

if r is located inside object 2, whereas ©iMiOj s (r, r) 
is only nonzero if r is located inside object 1. Thus, 
the integral over V2 in Eq. (|5T)|) can be extended over 
all space (without changing the result), if we restrict to 
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OiMiOj At, r). Then, the integral turns into a trace of 
this operator, and we finally have the exact result 



Rewriting 



in terms of the T operators, wc have 



HI 



(2) 



oo 

tJ 



dcu- 



Im Tr 



g k B T l — 1 



H 



(2) = 

2h 

7T 



-2ft 



du- 



ll e fc B T i - 1 



ImTrUl + GoTa) 



1 



du^— Tr { [Im[T 2 ] - T* Im[G ]T 2 

e k B Tl _ I 



1 - G Tl<KoT 2 
1 



1 - G0T1G0T: 



Jo [Im[Ti] - Ti Im[G ]Tl] G^ — 



■G [Im[Ti] - Ti Im[G ]TI] G^ 



T*G*TIG* 2 J' 
(55) 



1 - T^GpT^Go 



(56) 



The trace makes no reference to a specific basis, and 
the transfer is completely determined by the scattering 
properties of the two objects and the free Green's 
functions. This form may also be useful for cases where 
the T-operator is not known explicitly and has to be 
computed numerically pol |68| , as such methods are most 
powerful if reference to a specific basis can be avoided. 
As described below Eq. (|3"T[) , the trace in Eq. ([51)]) can be 
transformed into a trace over partial waves, in analogy 
to the procedure for equilibrium Casimir forces [7J. 
However, the resulting traces will have restrictions with 
respect to propagating or evanescent waves [see, e.g. 
Eq. (|102[) ]. In the following subsections, we shall prove 
the symmetry as well as the positivity of Eq. (|5GI) . 



where the radiation operator 
and 

W Q(3 = ^ 



is defined in Eq. ^1 



(58) 



u l-G T Q G Tiy 

The required symmetry is now apparent from Eq. (|57[) . 
as the trace allows a cyclic permutation of the operators, 
and furthermore we can take the complex conjugate of 
the expression since it is real. We have thus shown that 



H[ 2) (T) = H y 2>(T), 



(59) 



indeed holds for arbitrary objects. This allows to write 
the total heat transferred from object 1 to object 2, 



H 



l->2 



H 



(2), 



1 mj-i^m 

(e.g., the first) function, as [26 



in terms of one 



H 1 ' 



H?\T{)-H["(T 2 ) 



(2), 



(60) 



B. Symmetry of transfer 

The result in Eq. (|56|) can be used to prove the symme- 

(2) 

try of heat transfer: It is intuitively clear that H[ '(T) = 

H^p (T) has to hold since at equal temperatures the ob- 
jects should not exchange energy. While this is com- 
monly accepted, and has been shown numerically for the 
case of two spheres [37j , there is to our best knowledge 
no fundamental principle guaranteeing its validity. De- 
tailed balance cannot be invoked as even for Ti = T 2 the 
system is out of equilibrium if the environment is at a 
different temperature T env . For the case of two parallel 
plates, this symmetry is apparent from the formula for 
heat transfer [TJ, |36[ (in this case, it does follow from 
detailed balance as T onv plays no role). In order to prove 
the symmetry, we rewrite Eq. (|56[) as 

2ft f°° to 

h[ 2) = — duj^^ Tr{M*Wi 2 RiW^}, (57) 

n Jo e k B T i - 1 



C. Positivity of transfer 

We have shown in Section IIII CI that R a is a positive 
semi-definite operator, and the same holds for K* . With 
the property W 21 = Wj 2 for the operator in Eq. ([5"g[) . wc 
can write 

H[ 2) = —[ duj—^ Tr(lR;Wi 2 MiWl 2 } . (61) 

T Jo e^^" - 1 L J 

It is clear that Wi 2 RiW{ 2 is positive semi-definite, and 
thus also R 2 Wi 2 MiW{ 2 , because it is the product of two 

semi-definite operators. Thus the integrand in Eq. ([6T[) 

(2) 

is non-negative for any lj and the heat transfer H\ is 
non-negative, 

H[ 2) > . (62) 

This important proof has also to our knowledge not been 
presented before. As naturally expected, it shows that 
energy is always transferred from the warmer object to 
the colder one. 
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D. "Self emission and absorption, and the 
influence of other objects and environment 

Another important quantity is the heat emitted by ob- 
ject 1 in the proximity of object 2, which is for example 



relevant to the cooling rate of object 1. It is given by 
Eq. (f5U|) . with the integral taken over V\, and hence de- 



termined by the trace of < 



in Eq. (f5"Tj) , as 



J 



H 



(i) 



-2ft 



did z 

o e k B T i 



1 



ImTr (1 + G T 2 ) 



1 



1 - G TiG T 2 
I 



l [Im[Ti] - Ti Ini[G ]TJ] 



1 - G*T* 2 G*Tl 



(63) 



Note that — [as H a in Eq. ([37]) ] is positive if the ob- 

( 2) _____ 

ject emits energy. As it was the case for H\ ' in Eq. ([52"]) , 

the sign of is fixed. Interestingly, we can rewrite 
Eq. (|63p similarly to Eq. ([57)1 . where now the Green's 
function of object 1 appears, 



n Jo e k B T i 



-Tr{lm[Gi]Wi 2 MiWl 2 } . 

(64) 



As Im[Gi] is a positive semi-definite operator, we have 
proven that 



< 



(65) 



(2) 

holds for arbitrary objects. Having derived H± ' and 

-ff-P , we can now write down the total heat absorbed 
by object 2 for arbitrary temperatures T\, T 2 and T env . 
It is given by 

H { - 2 \Ti 1 T 2 , T cnv ) = H?\t x ) + Hf\T 2 ) + HW(T env ), 

(66) 

(2) 

where He n v(T cnv ) (which we do not give explicitly) is the 
radiation of the environment absorbed by object 2. Using 
Eq. (|17l) . we can express H^> solely in terms of and 



H. 



(2) 



H^(T 1 ,T 2 ,T cnv )= J2 Hi 2) (T a )-H^(T cl 



(67) 



a=l,2 



We stress that, using Eq. ([67]). the functions H& are 
sufficient to describe any heat balance for two objects, 
including the temperature of the environment. 



E. Generalization to N objects 

The generalization to N > 2 objects is 
straightforward- assuming that the composite T- 
operator of a collection of objects is known. We recall 
that the correlator C^f in Eq. (|4"!)|) is the radiation of 
object a scattered at all objects. Also for N > 2, this 



correlator carries Gq on its rightmost position; after the 
application of its inverse Gq 1 * in Eq. (|50[) . the operator 
on the rightmost position will be the T-opcrator of one 
of the objects [87J. Again, this is precisely because the 
final field correlator is always written in terms of the 
total currents on the objects, and expressing the corre- 
lation in terms of scattering operators allows to identify 
the sources with individual objects. In other words, 
E ■ (Gq U E*) in Eq. (J5DJI 

can always be decomposed as 







(68) 



where is the part which contains T% on the rightmost 
position, and is also the only term that contributes to the 
integral over the volume Vp . As before, the range of inte- 
gration can now be extended to all space since F„ (r,r) 
is only nonzero within Vp. Thus, also for N > 2 objects, 
we can write the heat absorbed by object j3 due to the 
sources in object a as a trace, 



* Jo 



duj — = — 

o fc H Tq 



ImTrF^. (69) 



The total heat absorbed by object /3 is then given by a 
sum over a, analogously to Eq. (1661) . 



H^({T a },T env ) = X)ff<f>(T a ) + #!«(T env ) , (70) 



where again we do not need to specify Henl(T env ) because 
we use Eq. (fT7|) to get 



H^({T a },T cnv ) = £ (HV\T a ) - H^(T cnv ) 



(71) 

We emphasize again that F^ can be expressed in terms 
of {T Q }, and Eq. is free of references to any specific 
basis. 
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VI. TRACE FORMULA FOR 
NON-EQUILIBRIUM FORCE 

A. Two objects 

Let us again start with two objects 1 and 2 at tempera- 
tures T\ and T 2 , respectively, in an environment at T env . 
We have shown in Ref. 271 that the forces on the two 



objects are not equal and opposite in non-equilibrium, 
and have to be derived separately. The force on one of 
the objects (say 2) is derived in close analogy to the heat 
transfer in Section [Vj This force can be found from the 
surface normal component of the Maxwell stress tensor 
a, integrated over the surface of object 2 (66j . 



F (2) = Re <b a n 
Js 2 



(72) 



with 



0-,,-lr 



where all fields arc evaluated at r. As for the radiation 
of a single object in Eq. (|37|) and the heat transfer in 



Eq. (|56j) , we can also derive a trace formula for the force 
by rewriting the surface integral in Eq. (f72j) as a vol- 
ume integral. Physically, the volume integral describes 
the Lorentz force acting on the fluctuating charges and 
currents inside the object. This leads by a straightfor- 
ward calculation to the following expression for the j th 
component of the force, 



j F 



(2) _ 



C 

4~ 



(73) 



V-2 



The force on object 2 has contributions due to all sources 
in the system [27], [29j]. We first consider which is 
due to the sources in object 1. It is determined by the 
correlator Cf c . We note that the splitting of the operator 
O1R1©} in Eq. (51) is helpful here too. Here only 
contributes to the integral in Eq. (|73l) in which case the 
range of integration can be extended over all space. We 
hence find for the force on object 2 due to the sources in 
object 1 



,(2) 



2ft f°° , 
— / dijj 

71" Jo 



e k B T t _ I 



ReTW V(1 + G T 2 ) 



1 



1 - G TiG T 2 



G [Im[Ti] - Ti Im[G ]T*] G^ 



1 



1 - T*G5T*GJS 



-To 



(74) 



We refer to this force in the following as interaction force [27j (this is the reason for the subscript i in Oj i ). The 

force on object 1 due to the sources in object 1, i.e., F^\ is referred to as the self force. F 2 2 ^ is found by exchanging 
indices 1 and 2 in the equation below. It is given by the self-part in Eq. (|52[) . and we have 



■in 



duj- 



e k B T 1 _ 1 



ReTr^ V(l + G T 2 ) 



1 



1 - G TiG T 2 



Go pm[Ti] - Ti Im[G ]TI] — 



(75) 



r 



The total force on object 2 is then given by a sum over 
all contributing (thermal and quantum) sources, 



£ [Fi 2 >(T Q )-F( 2 )(T onv )]. (77) 



F^({T a },T cm ) = J2 F^(T a ) + F^(T cnv ) + F { - 



(2) 



(76) 
(2) 

Fonv(Tenv) is the force on object 2 due to thermal fluctua- 

tions of the environment and F is the contribution from 

zero point fluctuations, i.e., the usual zero-tempcrature 

Casimir force. Using Eq. (|17l) the total force can be ex- 
(2) 

pressed in terms of {F a (T a )} and the equilibrium force 
F( 2 ' e «) as il 



F^({T a },T cm ) = F^\T cnv ) 



Equation (|77| has the advantage over Eq. ([76| that the 

(2) 

evaluation of the term Fo nv is not necessary (although 
possible with slightly more effort). Equations (|74)) . (f75|) 
and (|77p allow the computation of the non-equilibrium 
force between arbitrary objects in an arbitrary basis, and 
constitute another of our main results. 



B. Generalization to N objects 

The generalization to TV > 2 objects follows by con- 
sidering the composite T operator of the collection of 



objects. Using the decomposition of 
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Eq. (f68|) . we note that only contributes to the force 
on object /?, because it has on its rightmost position. 
After extending the range of integration to all space, we 
find 



F<?>(T )=** 

IT 



djjj- 



1 



ReTr 



1 



(78) 



The total force on object /? is given by Eq. ((77)) with the 
upper index 2 replaced by (3 and the sum running over 
all objects in the system. 



VII. PARTIAL WAVE REPRESENTATION 

A. Partial wave expansions of the free Green's 
function and the T operator 

1. Free Green's function 

In Sections IIII1 fVl and I VII we derived trace formulae 
for heat radiation, transfer and non-equilibrium forces, in 
terms of operators involving Go and T a . These formulae 



hold for any geometry. In this section, we present the 
derivation of the corresponding formulae in partial wave 
bases. As is the case for equilibrium Casimir forces (see, 
e.g., Refs. @, HI]), the traces of operators in the previ- 
ous sections will turn into sums over matrix elements 
with respect to partial wave indices, providing simple 
closed form equations for specific geometries. The ma- 
trix expressions below will have the general restrictions 
discussed in Refs. 0, H^. For example, two objects can 
only be described in a spherical basis if their enclosing 
spheres do not overlap, and a plane waves basis can only 
be used for two objects if they can be separated by a 
plane, see Fig. [T] 

In Ref. Q, the relevant functions (Green's functions, 
translation matrices and scattering amplitudes), were 
given for imaginary frequencies, as needed for equilibrium 
computations. In the present case of non-equilibrium, we 
have to evaluate them for real frequencies, which leads to 
some differences in the definitions as outlined in Appen- 
dices [5] to [EJ In contrast to Ref. Q , we choose to write 
the free Green's function without complex conjugations 
of waves, which yields an expansion that is manifestly 
analytic in the upper complex frequency plane, 



J 



<tO) 



^r)®E r ;s } (r') 

out ^j./^ 



(r) ® E° u ' 



if &(r)>£(r') 
if £i(r)<£i(r') 



(79) 



r 



Here [i runs over polarizations (electric and magnetic) as 
well as indices of vector-functions. £j is the 'radial' coor- 
dinate which gives rise to the two different pieces of 
the expansion in Eq. (|79|) . E res denote waves which are 
regular at the origin, and E ou * denote outgoing waves 
that are typically singular at the origin. The function 
<j{\i) is a permutation among the indices, which fulfills 
a(a(n)) = [i; along the lines of Ref. [701 , where, for 
example, in spherical waves o~({l,m}) = {l,—m}, i.e., 
the permutation changes the multipole index m to — m. 
See Appendices [Bl [Cl and [Dl for the expansions in plane, 
spherical and cylindrical waves. 



2. Definition of matrix elements of the T-operator and 
their relation to the scattering amplitude 

We define the matrix elements in accordance with the 
definition of the T operator in Eq. (|T9)) . The homoge- 
neous solution of the wave equation is the regular wave 
E™, s (r), and the scattering solution E sc is then 



E s ;,(r) = (l 



J T)E r ;, g 



E^(r) 



(81) 



Another important quantity that appears in the non- 
equilibrium formulae is the imaginary part of Gq. It is 
regular everywhere in space and hence can be expanded 
in regular waves 



where the second term is the scattered field. From this 
equation the matrix elements are obtained as [HI 



7^' — i 



d 8 r|d 3 r'E^ ) (r)T(ry)E7(r'). 



(82) 



Im[G (r,r')]=^E^(r) 



E 



(r'). (80) 



Here, the sum runs only over propagating waves. 



Note that the integrand involves no complex conjuga- 
tions (in contrast to e.g. Ref. Q). This ensures mani- 
fest analyticity of the matrix element 7^^' in the upper 
complex frequency plane. However, in most cases, the 
elements defined in Eq. (|8"2"j) are identical to those de- 
fined previously (e.g., for spheres (27| or for cylinders in 
Ref. (Hi). [H 
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FIG. 1: The results derived in the previous sections (Sees. IIIII - IVI|I are completely general, and apply e.g. to the configuration 
on the left hand side. The right hand side shows a configuration that allows representation in partial waves as derived in this 
section, because the enclosing spheres or ellipsoids do not overlap. 



The symmetry of T implies that the matrix elements 
7^' satisfy the condition 



Tim 7 — 7(7 



(m'Mm) 



(83) 



The S matrix [which we employ in Eq. (|97[) below] is 
defined as 



T 



S-I 



(84) 



3. Properties of partial waves 



The formulae for heat radiation, transfer and forces in 
Sees. IIIII IVl and IVII involve complex conjugates or imag- 
inary parts of the operators T Q and Go- This is in con- 
trast to the equilibrium force formula, see App. [X] and 
Ref. 0, which involves no complex conjugation or imag- 
inary parts. In order to relate, e.g., Im[T] to the matrix 
elements 7^' , we have to know the behavior of partial 
waves under complex conjugation. We note that the ex- 
pansion of Eq. (JTHJ) is not unique, as are the properties 
below. However, for any basis, the Green's function can 
be written in the form of Eq. (|79[) and the partial waves 
can be assumed to fulfill the properties below. Under 
complex conjugation we have for propagating modes 



E r ; s *(r) M epr. 



(85) 



For evanescent waves, complex conjugation involves a 
phase e l< ^ [H, 



E 



out 
<r(p) 



(r) = e*"E«s>) ^ 



E° 7 !U(r) 



. e -#ME° ut *(r) /xeev. 



(86) 
(87) 



The phase e^ f ' = e l< ^(f> can be easily found for any 
specific basis, as in Apps. [B] [C] and [D] Equations ([55)1 
and (|86p are identical for propagating modes where 
= i Using Eqs. CU), (HO]), d82]) as well as the rela- 
tions (|85p and (|86p . we can in a straightforward manner 
evaluate the operator expressions for radiation, transfer 
and forces, as demonstrated in the following subsections. 



B. Heat radiation 

In this subsection we omit for brevity the index a keep- 
ing in mind that all quantities refer to object a. 

1. Field correlator 

The partial wave representation of the field correlator 
C(T) for a single object in a cold environment in Eq. (|3"Tj) 
is now easily found and conveniently expressed in terms 
of outgoing waves, as 



Jo(T* — T) G* =iJ2 



x E° ut 



/vrt 
mm' 



■gout* 



"77 



and also 
G TIm[G ]T*GS 



E E W 



-pi out 



TTlOUt* 



(88) 



(89) 



Note that fj, and /jf run over all waves because the (near-) 
field correlator also contains evanescent waves. However, 
the index p in Eq. (|89p is restricted to propagating waves. 
The origin of this restriction is mathematically due to the 
restriction in Eq. (|8T)|) , and physically due to the fact that 
the environment radiation that enters the derivation of 
C via Eq. (f2"5)l . contains only propagating waves. It is 
useful to introduce the projector on propagating waves, 

n^ = <W<W- (so) 

With this definition the net correlator is 



C(T) 



— 1$ 



x E° ut <g> E? ut * 



(91) 

We also define the matrix version of the radiation oper- 
ator K in Eq. ([3Tj). 



— i$ 



(92) 
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which yields for the emitted field 



Heat transfer 



out* 



(93) 



C(T) = a(T)^7^E° ut ®E 

n,n> 

2. Emitted energy 



The operator-trace for the emitted energy in Eq. (|37[) 
is readily written as a trace of an (infinite) matrix. For 
example, the first term in Eq. (|37j) reads 



Tr{Im[G ] Im[T]} = 
/ / d 3 rd 3 r' ]T E^ ) (r')Im[T(r',r)]E^(r) 

= - ^ Re T MM = -Tr pr Rc T, (94) 



where we have also used the symmetry of T. In the 
last equality, we have defined the trace over propagating 
waves of the matrix 77m' ■ The second term in Eq. ([371) is 



treated analogously, and we obtain for the emitted energy 
in a partial wave basis 



H = — — f Tr pr {Rc[T] + TT f } 

K J e k B T — 1 



(95) 



In the second term, we defined Tr pr [TT^] = X^/x' \T~w' | 2 > 
where both fi and fx' run over propagating modes. 

Given our definition of the matrix form of the radiation 
operator in Eq. (|92p , we can also write the emitted energy 
as 



H 



— / 

7T 



- 1 



-Tr pr TZ. 



(96) 



There is an alternative representation of the emitted en- 
ergy in terms of the scattering matrix S known in liter- 
ature [62|, HH , which also follows from Eq. ([84)) as 



2^ J pk R T . 



1 



-Tr pr [1 - SS^] 



(97) 



In Ref . [62| , Eq. ([97} was derived by an entirely different 
route, starting directly from partial waves expansions and 
using identities for vector waves. 



(2) 

Let us transform the heat transfer expression H\ in 
Eq. (|56p into a partial waves basis. First we note that 
Eq. (|56j) contains Go in two different ways, (i) Im[Go] 
connecting T operators of the same object. While the 
free Green's function is singular at the origin, its imagi- 
nary part is regular and can be expanded using Eq. (|80l) . 
This case where Go is sandwiched by T operators of the 
same object does not occur in the equilibrium Casimir 
formula [7f. (ii) The remaining Go in Eq. (|56[) are similar 
to equilibrium as they connect Ti and T 2 . To expand 
the latter we use the techniques presented in detail in 
Ref. , expanding the outgoing waves in Eq. ([75} in the 
coordinate system of the other object as 



TTlOUt 



E 



Wg(X a p)E7(r a ), 



(98) 



where r Q and denote the same position measured rel- 
ative to the origin of systems a and ft respectively, and 
Xq,^ = v a — r« is the vector connecting the two coor- 
dinate origins @, HH, IzH- We can then write the free 
Green's function as 



)(r 2 ,ri) 



■ ^E7(r 2 )®E^ 



(99) 



Here, and r 2 are measured in the coordinate system of 
objects 1 and 2 respectively, and are located on the cor- 
responding object. Applying the expansions in Eqs. ([80)) 
and (|M[) to Eq. ([55|) . we arrive at the expression for the 
heat transfer rate in matrix form. 



1. Spherical Basis 

We start with the spherical basis, as it is most useful 
for studies of compact objects. It also allows for the most 
concise representation as it docs not contain evanescent 
modes. The expression for heat transfer is 



(2) 



2// 



dw- 



-Tr 



s T i - 1 



r; + r 2 



■7?T 2 



1 



2 - UT1UT2 



U 



Tl + Ti 



wH, (100) 



where the trace involves summing over all indices 
{I, m, P}, of spherical waves; and the adjoint of U is to be 



taken with respect to all these indices, i.e. = U 1 ^*, . 
For example, one of the terms (linear in the T matrices) 
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reads explicitly Tr^UTi^] = V^Uf^Tx^^Up,, 
with sums over all indices. Equation (|100[) is valid in 
spherical basis, but the objects described can be of any 
shape. For homogeneous spheres, the matrix Tp^im is 
diagonal in Z, m and P, and Eq. ()100[) simplifies further 
(see Sec. lXAj) . 

2. Arbitrary Basis 

In general, the wave expansion in Eq. (|79[) contains also 
evanescent waves, which behave differently under com- 



plex conjugation, sec Eq. (|87j) . The consequence is that 
Eq. (jlOOp is modified and contains factors of e 1 ^ in some 
places [compare to Eq. ([91]) ]. The result does however 
take a simple form in terms of the redefined matrices 

r^= e -^T^, v^,=U^,e 1 ^'. (101) 

Because Im[Go] contains only propagating waves, the 
projection ±l pr of Eq. ([90]) appears between the matri- 
ces of the same object, such that the heat transfer in an 
arbitrary basis is 



H 



(2) 



2h 

7T 



dw- 



-Tr 



e k B T± _ I 



To 



VT\VTl 



(102) 



We note that the matrix r is in general non-analytic in 
the upper complex frequency plane, whereas T is man- 
ifestly analytic. Equation (|102[) includes Eq. (|100[) as a 
special case since in the spherical basis r = T, v = hi 
and IP r = I. 



D. Non-equilibrium Force 

The traces in Eqs. (|74|) and (|75|) can also be trans- 
formed to sums over partial waves. The discussion be- 
fore Eq. (|100p holds in close analogy here, too. Note that 
Eqs. (|74[) and (|75|) are directly for the force, whereas 
in equilibrium one starts with a free energy Q. In 
App. [A] we demonstrate that the force obtained in our 
non-equilibrium formalism, if applied to equilibrium, can 
be integrated to yield the free energy; while more gen- 
erally Eqs. (|74"|) and (jT5j) cannot be thus integrated. In 
order to evaluate the gradients in Eqs. (fT4"|) and ([75)1 . we 
note that they act in two different ways on Go- In the 
first case, the gradient acts on a Go that connects Ti 
and T2. After using Eq. (|99|) . this action turns into the 
derivative of U(d), 

V r2 G (r 2 ,r!) = 

= -*E V ^^( X 2i)E7(r 2 )®E^ ) (r 1 ), 

= *E(p^ 21 W ( X ^) E 7M ® E^j(n). (103) 
I 



In the last expression, we have introduced the infinitesi- 
mal translation operator 

P w ' = -V a V w ,(a)| a=0 , (104) 
where V,,,,, describes the translation of regular waves in 
the same basis by 

Er ; s (n) = E w^ 1 ^ + a )- ( 105 ) 

The two representations of Eq. (|103j) can be equivalently 
used; we present formulae below in terms of the second 
representation which allows for a more compact notation. 
In the second case, the gradient acts on the free Green's 
function connecting T2 with itself (a combination which 
does not appear in calculations of equilibrium forces). 
This gradient can also be easily expressed in terms of p, 
[where in contrast to Eq. (|103j) . both points are measured 
in the same coordinate system] as 

V r Im[G ](r,r')= E E P^^V « ® ^ 

(106) 

Apart from these derivatives, the result is in close analogy 
to the heat transfer in Eq. (|102|) . Here, for brevity, we 
give only the result for arbitrary basis, see Eq. (|10ip , [In 
particular, in the spherical basis, t = T ", v = U and 
n pr = X hold, further simplifying the expression] 



r (2) 



2h 

7T 



duj- 



Im Tr 



g k B T l — 1 



T 2 t p + T 2 t np r P np r T 2 



VT\VTi 



(107) 



15 



The self contribution of Eq. (|75| is also readily written in a partial wave basis as 



2h 



duj- 



-,k B T 1 



Im Tr <^ [p v t 2 v + II pr pII pr ] 



I — vW^vWl 



(108) 



For the meaning of the matrix multiplications in 
Eqs. (fTU7|) and (fTT75|) . sec the discussion below Eq. (jTUOT) . 
For, e.g., homogeneous spheres, the expressions simplify 
further as described in Sec. IXI Al 



VIII. KEY DIFFERENCES BETWEEN 
EQUILIBRIUM AND NON-EQUILIBRIUM 

In contrast to equilibrium force calculations, the ex- 
pressions for heat transfer and non-equilibrium interac- 
tions are non-analytic in the upper complex frequency 
plane (e.g. due to the presence of adjoint quantities) and 
have to be evaluated for real frequencies. This restriction 
reflects the fact that non-equilibrium quantities can be 
strongly influenced by small changes in the resonances 
of the dielectric functions [U Hfl, effects of which are 
marginal along the imaginary frequency axis. While non- 
equilibrium effects are hence richer in their phenomenol- 
ogy, they are also harder to evaluate numerically due to 
oscillatory behavior of the functions involved. 

Setting aside the issue of convergence of the series in 
Eq. P5]) at close proximity, the evaluation of heat transfer 
and interactions can be simplified in the limit where the 
separation d is much larger than the size of the objects, 
where the following one reflection approximation 

O««l + G T a , (109) 

becomes asymptotically exact. This is because most of 
the waves involved will not scatter twice at the same ob- 
ject. We will use this approximation in Sees. IXl and IXj 
to derive analytic expressions for the cases of a sphere in 
front a plate, and for two spheres. In these cases we ob- 
serve a difference in convergence of multipole expansions 
between equilibrium and non-equilibrium situations. For 
equilibrium Casimir interactions (see e.g. Refs. {EI l72j). 
the convergence of the multipole expansion is governed 
by the ratio R/d s , because the dominant wavelength is of 
the order of the surface to surface separation d s . For heat 
transfer [37} and non-equilibrium forces [27|, the conver- 
gence is governed by both R/d s and R/Xt, and the max- 
imum multipole order to be included is roughly given 
by the larger of the two numbers (see Refs. [33, lZ3[ f° r 
more detailed discussions regarding heat transfer) . Even 
if R/d s is small, one might still need many multipoles 
if R/Xt is not small. For R <C At, however, there is a 
well defined asymptotic large d expansion, equivalent to 
the Casimir Polder limit plus higher order corrections, in 
equilibrium For example, Eqs. (|135[) or (|147|) (and 
in general all our equations for R <C Ay), are asymptotic 



large d expansions, in the sense that the given orders in 
inverse d have no corrections from higher order reflec- 
tions. 



IX. EXAMPLES FOR HEAT RADIATION 
A. Radiation of a plate 

The radiation of a plate, i.e., a semi- infinite planar 
body occupying the space z < 0, has been extensively 
studied by many authors (l3l . |64 |. The main empha- 
sis is on the radiated energy, whereas the correlator of 
the emitted field, as discussed in Sec. Mil can be found, 
e.g., in Ref. [l4j]. The simplicity of the plate geometry 
(labeled by p) enables easy derivation of formulae, e.g. 
from Eqs. ((57)) and (|3T)]) . We will use this geometry to 
demonstrate two things: First, we show that Eqs. (|9"Tj) 
and (|95|) can indeed be applied to planar geometries, 
which is slightly less obvious than for spheres as below. 
Second, we show that it is not strictly necessary to use a 
basis that satisfies the conditions of Eqs. (j85l) (|86l) . It is 
sufficient to have a transformation from the desired basis 
to a basis that docs fulfill Eqs. (|g5 j) - ([8l))) . 

We first consider the emission of a plane-parallel di- 
electric slab of finite-thickness in the region — I < z < 0. 
We compute the correlator of the electric field for any 
two points r and r' outside the slab using the general 
formula of Eq. (j9"Tj) . For simplicity, we consider the case 
when both r and r' are on the side z > 0, so that the 
field can be expressed in terms of right traveling waves 
E< R U p,ki ■ As described in App. [5J Eq. (jHI]) is not appli- 
cable to such waves, but holds for the waves of definite 
parity, E° u p k ^. However, we can still use Eq. (|9T|). not- 
ing that the two sets of waves are related by the unitary 
transformation 

E° u p,k,= E E $,k <V (HO) 

j=L.R 

This transformation follows immediately from Eq. (|F38|) . 

Oj S = -j=[5 s ^ + 8j, R + <5 s .+<5j,l + iS s -Sj t R - iSg-Sj^)- 

(111) 

Because E° : u p ki (r) = for z > [see Eq. (|B5)) ]. we have 




x Kr,p,^,R,p'M' ± E£f P , k± (r) ® E^, )ki (r') , (112) 
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where we have introduced the transformed radiation op- 
erator 

K = OTZO ] = -- (f e _< * + e*+f f ) - fU^f 1 , (113) 

given in terms of the transformed T-operator, T = 
OTO\ and = Oe^O^. In contrast to e^, is 
not diagonal, and we find by explicit computation that 

x[e pr %,-i(l-<5^)e ev ] , (114) 
where we have used the step functions 



e pr = Q(u/c-k±), 



G(kx-u/c). (115) 



vanishes and rp® approaches r p given in Eq. (|B12[) . The 
emission to the right is then 



El = A / 

A 2tt 



e k B T p — 1 J k±<^ 



d 2 k 



(27T) 



(120) 



A slab made of a perfectly reflecting material (|e| — ^ oo) 
does not emit energy as in this case the Fresnel coeffi- 
cients approach unity (and the transmission vanishes). 
This is a manifestation of Kirchhoff's law: since a plate 
with perfect reflectivity does not absorb electromagnetic 
waves, it can also not radiate. The blackbody limit is 
obtained by letting r p — > in Eq. (|120[) . in which case 
Hp approaches the Stefan-Boltzmann law [HI, [HJ . 



According to Eqs. (|8T|) and (|82[). the matrix T is re- 
lated to the Fresnel reflection and transmission coeffi- 
cients, rp and tp for outgoing waves to the right of 
the slab, by 

7R,P,k±,R,P',k' ± 



<W,(27r) 2 <5( 2 >(k ± -k^) ^ 



(R) 



Tr 



R,P,k±,L,P' ,k' 



5p P ,(2ir) 2 5 (2 Hkx-k' ± 



„(«) 



(116) 
(117) 



Both t^ fl) and r { p ] (see e.g. Ref. [zl p.299]) depend on 
the thickness I of the slab. Substitution of the above 



matrix elements in Eq. (|112j) gives 
d 2 k A 



C p( T p) = a(T p ) 



p 



(27T) 2 



i 



.(«)|2 



i4 fl) i 2 



e pr 



E 



P,kj_( r ')- 

(118) 



The first term in Eq. (|1 18[) describes propagating waves 
which carry energy emitted by the slab. The second term 
corresponds to evanescent waves, which do not contribute 
to the energy emitted. In the situation of heat trans- 
fer between multiple objects, it is the evanescent waves 
which lead to a strong increase of transfer at close sepa- 
rations [Hj]. Equation (pj)) for the total emitted energy 
per surface area A is also readily evaluated to give 



B. Radiation of a sphere 

Another known result is for heat emitted by a 
sphere [5lT |. which we re-derive here in our notation 
for completeness. (We also provide the field correlator, 
which seems not available in the literature.) The spher- 
ical basis (see Appendix [C] for details) is naturally the 
most appropriate, and Eqs. (j9"Tj) and (|95[) are readily eval- 
uated. For simplicity, we consider a homogeneous sphere 
of radius R. The matrix 7^' is diagonal [see Eq. (|C5[) ], 

C s (T) = -a(T) J2 [Re7I P + |7;T] 

P,l,m 

Ept n (r,6,cj>)®Ep a Z(r',9' >( t>'). (121) 

The trace over the index fi now involves sums over P, 
I, and m. The rate of energy emission by the sphere is 
found from Eqs. (|3"?| or (f9"5")) . as 

H s = -B f i;-^ £ [RcT, p + \T t p \ 2 ] ■ 
ir Jo c*bt -l Ptl 



(122) 



Equation (|122[) can also be obtained from Eq. (|121[) by 
integrating the Poynting vector over the surface of the 
sphere E s , where the following identity is useful, 



Im 



(E^OR, M) x V x E°pf*(R, 0, 0)) • f = 1. 



H p (T p ) 



dio- 



27T ,/q gfc R 5 



d 2 k ± 



(119) 



In view of Eq. (|1 18[) , H p represents the sum of the emis- 
sions of the two faces of the slab, the emission of the right 
(left) face being given by the first (second) pair of round 
brackets. For an infinitely thick plate, the transmission 



In the limit of perfect conductivity (or reflectivity), the 
emission of the sphere vanishes because one has 



lim Re[T t P ] = - lim \T t 

\e\— >oo \e\— too 



P\2 



(123) 



Equation (|123[) is a general property of scattering oper- 
ators, which is connected to the fact that the 5-matrix 
in Eq. (|84[) becomes unitary in this limit. Objects with 
unitary £ do not absorb energy and hence cannot radiate 
heat [EllillHlzl. 
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FIG. 2: Energy emitted by SiC>2 (upper curves) and Gold 
(lower curves) spheres at T — 300K in a cold environment as 
a function of radius R, normalized by the Stefan-Boltzmann 
result (Eq. (|124p with e(T) = 1). Solid lines correspond to the 
exact result from Eq. (|122|l (for Si02, the data are from [2(|); 
different approximations valid for small spheres are depicted 
by: Long dashes correspond to Eq. (|127|1 . short dashes to 
Eq. ffT29]l omitting the term |7I P | 2 , and dots to Eq. ([T29J. 



Figure [5] shows the radiation of a sphere for two differ- 
ent materials, SiC>2 and gold, as examples for dielectrics 
and conductors, respectively. The SiC>2 sphere (where we 
have used optical data) has been analyzed in Ref. [26| : 
If R is much larger than the thermal wavelength and the 
skin depth S = c/(1yh ^/efiuj), the emitted heat becomes 
proportional to the surface area of the sphere, and can 
be written as 



lim H s 



4irR 2 aT 4 e(T) 



A 



(124) 



with a = TT 2 kg/(60h 3 c 2 ). Here, e(T) is the emissivity of 
the plate, which can be derived from Eq. (|120p as shown 
by the last equality in Eq. (|124j) (e(T) — > 1 corresponds 
to a black body). 

In the opposite limit, where R is the smallest scale, 
the emission is proportional to the volume of the sphere, 
and the normalized curve in Fig. [5] is linear in R. In this 
limit, one can use the expansion of the T-matrix for small 
R* = ^, as 



3(e + 2) 
4(e-l) 2 _^ 6 



2 . 2-3 £ + £ 2 (l + M ) E , 5 



9(2 + eY 



5(2 + e) 2 

o (r* 7 ) , 



(125) 



and accordingly for T^ 1 with e <-> /i. The first term in 
Eq. (|125[) is commonly attributed to the dipole polariz- 
ability 



e - 1 



K 6 . 



(126) 



For jU = 1, the radiation is then given by 

/ duj— ^ 

C°7T 



lim H s 

fl«{A T ,<5} 



Ima. 







1 



(127) 



While this is a good approximation for small SiC>2 
spheres, it can have a very limited range of validity for 
other materials, as can be seen from comparison to gold 
in Fig. f2J For the dielectric response of gold we have used 
the Drude model 



£Au(w) = 1 



uj{uj + iuj T ) ' 



(128) 



with uj p = 9.03 eV and uj t = 2.67 x 1CT 2 eV. For large 
R, the blackbody limit is approached, but with a much 
smaller emissivity compared to SiC>2 [t(T) vanishes as 
l/-v/e, as can be seen from expanding H p in Eq. (|120[) 
for large e]. However, the limit of Eq. (|127|) is not ap- 
proached for the physically accessible radii, as can be seen 
in Fig. [2] The relevant skin depth for gold is roughly 
20 nm, and the exact curve approaches the result of 
Eq. (|127[) for R rj 1 nm. However the radiation from a 
sphere of R <S^. At for all materials can be approximated 
by restricting the sum in Eq. (|122[) to I = 1, as 



lim H s 

-R«A T 



/ dio - 



1 p 



[Re 7/ 



IT/?] • 
(129) 

As indicated by the dotted lines in Fig. [21 this result holds 
for both dielectric and conducting spheres, accurately 
describing the emission up to roughly R/Xt ~ 10%. 
In Fig. [2j we also show the result after omitting the 
quadratic term in Eq. (|129j) . an approximation used in 
Ref. [13] in order to obtain short equations for the non- 
equilibrium interactions between spheres and for a sphere 
in front of a plate. While we expect this additional ap- 
proximation to have a smaller range of validity compared 
to Eq. (fT2"9l) . we find that it is as good as Eq. (fT27)) for 
Si0 2 , and much better than Eq. (fT27)l for gold. Addi- 
tionally, in contrast to Eq. f)12Tf) . it can be extended to 
magnetic materials. See Table U for a summary of the 
various approximations and their limits of validity. 

Lastly, we note that when normalized by the volume of 
the sphere, the radiation of SiC>2 decreases monotonically, 
while that of gold has a sharp maximum around R = 
100 nm, a feature which is not captured by the dipole 
approximation of Eq. (|127|) . 



C. Radiation of a cylinder 

An infinite cylinder is the last shape for which heat 
radiation can be computed analytically. The emitted 
energy was originally described in Ref. [5fjj . and more 
recently re-derived, and numerically analyzed in detail, 
in Refs. [Ill H3, [55[, paying special attention to the po- 
larized nature of the radiation. The correlation of the 
emitted fields was also recently studied in Ref. (30j . 
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-(21 + 1) y, D (Re77 p + |77 p | 2 ) « 


Range of validity 


-3Ep(Rer 1 p + |r 1 p | 2 ) 


-R < A T 


-3£ p Re7f 


l^/peli? < At 




\Ve\R < A T , m = 1 



TABLE I: Range of validity of different approximations for 
the radiation of a sphere. 



In cylindrical vector basis (see App.|D]for details), the 
matrix T^n/tor a cylinder is not diagonal in polariza- 
tions HI I!, [zH, but symmetric. With Eq. (JHU), the 
correlation of the emitted field is then readily evaluated 
as 



C C (T C ) = -a(T c ) E / ^ E 

n =-°°-oo 



E°4(r)®E* = (r'). 

(130) 



The quadratic term carries the step function O as it only 
spans propagating waves with — — to — , whereas the term 
linear in T contains both propagating and evanescent 
waves. The phase factor of the linear term is the real 
function 

e**- =[e(w/c-|A,|) +(-l)(" +1 )6(|fc 2 |- W / c )]. 

(131) 

The energy radiated by the cylinder (2(| [Ho|, HH follows 
from Eq. ([95]) [or Eq. $J$\, and reads 



r p n— — oo 



UJ / c 

j d ^{Mrm + Kl\ 2 + \r n p l\ 2 ), (132) 

—a; /c 

where P = M if P = N and vice versa. We note that 
Eq. (|132l) can also be found from Eq. (|130p by computing 
the Poynting vector, where the following relation is useful 

Im i / [E^ fc< (r) x V x E°J% k (r')] • p = <W . 

(133) 



In Rcf. [55|, the emission of a gold wire was found to be 
very large at R ~ 20 nm, where it exceeds the Stefan- 
Boltzmann value by a factor of 10. This feature is not 
found for the gold sphere in Fig. [51 which we attribute to 
the fact that the gold wire allows for unrestricted electric 
fields along the wire, but the sphere does not. 
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FIG. 3: The system of two spheres in thermal non- 
equilibrium. 



X. EXAMPLES FOR HEAT TRANSFER 



A. Two spheres 



The heat transfer between two spheres was studied nu- 
mericallyin great detail in Refs. [37l.[38|. and analytically 
in Rcf. [36j . Our Eq. f| L00[) provides a compact exact rep- 
resentation for this transfer in a spherical basis. Here, we 
present a more explicit evaluation for the case of homoge- 
neous spheres, where the matrices 7^' are diagonal [see 
Eq. (|C5p j. For simplicity, we consider only the one reflec- 
tion approximation. Assuming that the spheres are small 
compared to the thermal wavelengths, we then provide 
an asymptotic expansion valid at large separations. This 
extends the results in Rcf. [36} to spheres of arbitrary 
material. The spheres have radii Rj (j = 1,2), complex 
dielectric and magnetic permeabilities Sj and fj,j , at tem- 
peratures Tj, with their centers separated by a distance 
d, as in Fig. [3] The trace in Eq. (|100[) simplifies further 
since 7^' is diagonal, and the expression for transfer 
from sphere 1 to sphere 2 reads (in the one reflection 
appr oximat ion ) 



lim 7J 1 ^ 2 (T 1 ,T 2 ) 



2h 

7T 



£ k B T l — 1 £ k B T 2 — 1 



E (MT 1 p 1 ] + \t 1 p 1 \ 2 )(mt 2 



PP'W 



F" 
2.1' 



\uf x 



IP'Pml 



(134) 
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Equation (|134[) contains the translation matrices U for 
spherical waves (see App. IE 2\ . If the radii are small 
compared to the separation as well as the thermal wave- 
lengths, the above expression can be evaluated using 
Imax = 1 = I'max ( see a l so Fig- [2] and Table HJ. This 
corresponds to an asymptotic large d expansion of the 



transfer (denoting 7~ p 



ters. Equation (| 135|) requires only {d, At,} 3> Rj, and 
thus holds for all materials. Also the limit of perfect 
reflectivity (|e| — > oo) is only captured correctly when 
including the quadratic terms in Eq. (j!35[) . as only then 
the transfer asymptotically approaches zero [compare to 
Eq. dUSD] [9l|. 



lim H 1 ^ 2 



2h 



duj 



e k B T l — 1 C k B T 2 — 1 



x £ (Re[7f ] + |7f| 2 ) (Rc[T/'] + \Tf\*) 



p,p' 



9c 2 



9c 4 



27c 6 



2u 2 d 2 2w 4 d 4 2uj 6 d 6 



5p, 



(135) 



This equation has a common limit with Eq. (62) in 
Ref. [36j , if we restrict to the first term in Eq. (|125[) and 
linearize Eq. (|135[) in a. As we have shown in Fig. [2] and 
Table HI such a replacement is only valid if the spheres 
are non-magnetic and small compared to the skin-depth 
of the material; it is in general not valid for conductors 
for which the skin-depth is of the order of a few nanome- H pr + H ev as in Eq. (|118[) . We have 



B. Sphere and plate 

The heat transfer between a sphere and a plate was 
considered numerically in Refs. (2(| [3^, H(| and analyti- 
cally in Ref. [3(|. Here, we provide the result in a one- 
reflection approximation, including the large d expansion 
for small R/Xt- The exact transfer given in Eq. (|102[) 
is more complicated because every reflection involves an 
integration over wave-vectors. We consider the system 
shown in Fig. @] with the sphere described by R, e s and 
fi s , and the plate by e p and /x p . The sphere-center to 
surface separation is denoted by d. 

The heat transfer can be split into contributions from 



propagating (H pr ) and evanescent (H ev ) waves, H 1 




with pr and 6 ov given in Eq. (|1 15|) . If the sphere is small 
compared to the thermal wavelength, we can restrict to 
the terms with I = 1 (see Fig. [3] and Table H]), to get the 
asymptotic result for large d as 
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The distance independent term due to propagating waves 
now becomes (denoting T p = 7^), 
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Keeping only the terms linear in T p in Eq. (|138p . and 
simplifying these by restricting to electric and magnetic 
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FIG. 4: The system of a sphere in front a plate in thermal 
non-equilibrium. 



dipole polarizabilities, Eq. (9) of Ref. [26j is reproduced. 
The evanescent part can be analyzed in the two limits of 
d -^i Ay and d At- If d is the largest scale, i.e., for 
d ^> At ^> R, the evanescent part of the transfer decays 
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as d 2 (setting fi p = 1), 
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In the opposite limit, where the separation is much 
smaller than At but still much larger than R, i.e., for 
At q S> d^- R, we find a ~ d~ 3 fall-off, and 



3hc 

2nd 3 



~~2 







1 



1 



e f 



Im 



l + e % 



(Re[T N ] + \T 



N\2\ 



(140) 



Equation (|140j) is consistent with the first term in 
Eq. (28) of Ref. [H (after setting T 



JV 



3c- 



■a, and 

restricting to terms linear in the dipolc polarizability a). 
We note again that only the expressions in Eqs. (|138[) 
to (|140[) . including quadratic terms, are valid for small 
spheres (R <C At) made of any material. 



XI. EXAMPLES FOR NON-EQUILIBRIUM 
INTERACTIONS 



A. Two spheres 

Let us consider the Casimir interaction of the two 
spheres in Fig. [3] (again, with Rj, Ej and fXj) at temper- 
atures Tj in an environment at T cnv . The centers of the 
two sphere, 0\ and 2 , are separated by 0\ — 2 = dz. 
We derive the total force F^ 2 * 1 acting on sphere 2 (F^ is 
then found by interchanging indices 1 and 2 everywhere) . 
This force has three contributions in Eq. ([77]) : The equi- 
librium force for the two spheres evaluated at T onv , and 
two contributions due to the deviations of T\ from T env 
(Ff } ), and of T 2 from T cnv (F^ 2) ). Physically, these forces 
follow from the heat radiation of sphere 1 and 2, respec- 
tively [see Eq. ([T2"2"j)]. Equations (fTDTJ) and ([105)1 give 
the exact expressions for these forces between arbitrary 
objects in a basis independent form. For homogeneous 
spheres, these equations simplify due to the diagonality 
of the matrix [see Eq. ()C5[) ] , as well as the absence 
of evanescent modes. Here, we give a more explicit form 
of the result in a one reflection approximation, including 
all wave indices. The force points from center to center, 
and we define it to be positive when it is directed towards 
the other sphere (attractive). From Eq. (|107[) . we have 
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with P = N if P = M and vice versa. In the second 
expression, we have expressed the derivative of the matrix 
elements of U in terms of other elements of it, e.g., by 
using d d U = -p z U and d d V(dz)\ d=a = -p z , with p z 
given in Eq. ()E6)I . In doing so, the functions 
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1(1 + 1)' 



(143) 



l + l 



(2Z + 1)(2Z + 3) 



appear. The first representation, Eq. (|14ip . can be evalu- 
ated numerically or analytically. We show also the second 



representation, because it is found when performing the 
calculation as described in Ref. [27j, i.e., by evaluating 
the Maxwell stress tensor instead of using the trace for- 
mula of Eq. (|107j) . The product of two spherical waves 
gives a net force on the sphere only if they differ in polar- 
ization or by 1 in the multipole order / |76| , as a spheri- 
cally symmetric field generates no net force. 



For the self force, we evaluate Eq. (|108p in the one 
reflection approximation, also writing the matrices more 
explicitly, such that 
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In the second expression, we have again re-expressed 
the d-dcrivative in order to find the form which natu- 
rally arises when considering the stress tensor, as done in 
Ref. [13. We note that in both Eq. (fT43)l and Eq. fT35]) , 
none of the U matrices is conjugated, in accordance with 
Eq. (|108|) . This can lead to an oscillatory behavior of the 
self force as a function of d [27j ■ 

When the spheres are small compared to the thermal 
wavelength, R <C At, we can restrict the partial wave 
sum to the dipole moment (1 = 1) in all T matrices. 
This results in an asymptotic large d expansion which 
is the non-equilibrium counterpart of the Casimir Polder 
limit, as (denoting again Tf = Tf l=1 ) 
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This expression is identical to Eq. (6) in Ref. [27j , except 
that we have here included the terms quadratic in T . As 
shown in Fig. [5] and Table HI Eq. (|147|) holds for any ma- 
terial, while Eq. (6) in Ref. [27| holds only for materials 
with sufficiently small real and imaginary part of e — 1 
and \x — 1 . 

For large separations, F± decays as d~ 2 and is re- 
pulsive. This originates from momentum transfer to the 
second sphere via absorption or scattering of photons. 
The remaining terms in Eq. (|147|) . with higher powers 
in 1/d, are (in most cases) attractive. Similarly, the self 
force F. 

F^ 
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2 (T2) is expanded for {d, At 2 } 3> Rj, 

h 



as 



en 



dui- 



1 1 



(7f r p ) 



9c 2 



-1 p 



(Re[T/ 



2 d 2 



Tf . 144c 5 



j 5 d 5 



27c 3 



162c b 



(7? 



81c 7 



r 2 p | 2 ) 



T p 
— ) 



72c 4 

^ 4 d 4 



; 7 d 7 



(148) 



(This expression is, up to quadratic terms, identical to 
Eq. (7) in Ref. [27|.) We emphasize again that, in con- 

(2) 

trast to F{ , this term can oscillate as a function of d at 
a scale set by material resonances (see Ref. (27j)- 

We now repeat the low temperature expansion for di- 
electrics, which was given in Ref. [27|. The leading low 
temperature behavior of the force for insulators can be 
derived by requiring At 3> Ao , where Ao is the wavelength 
of the lowest resonance of the material. The dielectric 
functions and polarizabilities are then expanded as [66[ 
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(149) 



atj = a Q j + ia i0 j + 0(lj 2 ), (150) 

with £Q t j, Xi n j, cto.j an d ctioj = 3i? 3 /(eo.j +2) 2 real. For 
Atj ^> Ao, the interaction term is then given in closed 
form by 
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The self force F% does not oscillate to lowest order in 
temperature and takes a more complicated form. In the 
limit where d is the largest scale, we have 
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While in this range of d the force F 2 is independent of 
temperature, it vanishes as T 2 — > since with At 2 the 
largest scale (At 2 3> {d, -Rj, Ao}), one has 
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(152) 



lim Fr 



(2) 



d»Pj 2 d 7 A T 

n(2) 



which is identical to F^ in this limit, with indices 1 and 
2 interchanged. 



B. Sphere and plate 

The force on a small sphere in front of a plate, (see 
Fig. 0]), has been studied in Ref. [27| for various com- 
binations of temperatures (0 K or 300 K). The forces 
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on the sphere and the plate are not equal and opposite 
in non-equilibrium. We focus here on the force expe- 
rienced by the sphere, which is directed normal to the 
plate. The force is defined to be positive when the 
sphere, centered at z = d, is attracted to the plate (oc- 
cupying the space z < 0). The total force on the sphere 
given in Eq. (|77[1 contains three contributions: The equi- 
librium force at the temperature of the environment, and 
two contributions from the deviations of T s and T„ from 



( s ) is) 

T env , denoted by Fs and F p ' . These are most con- 
veniently expressed in terms of the formulae translating 
plane waves to spherical waves (see App. IE ip . 

The interaction force is naturally split into contribu- 

(s) 

tions from propagating and evanescent waves, Fp = 

(s) (s) 

Fpjr + FpJ v . In a one reflection approximation the for- 
mer is independent of separation, and 
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The unit-step functions pr and GV defined in Eq. (|115|) project onto propagating and evanescent modes, respectively. 
The self force due to the sources in the sphere comes from the radiation of the sphere at temperature T s . In the one 
reflection approximation, we find 
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where the function g(k±) contains projectors on propa- 
gating and evanescent modes, 



= (-i) m+i (i - 2<w)e pr + e c 



(156) 



Fs contains no separation independent term and van- 
ishes for d — > 00. Let us perform the simplifications that 
arise for the case R <§; At, where the restriction to dipole 
order, / max = 1, yields again the asymptotic large d ex- 
pansion. We obtain for the Casimir Polder limit of the 
interaction force 
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fpr — [ 



,/c 



^^(l-l^l 2 )^'] 

P,P' 



Re[T p T p *]5 PP , 



(158) 



lev — 2 1 



oj I c 



E ( Im V\ 



- 1 Im[T P ] - lm[T P T P *} 



lm[r p ] lm[T P ] 



(159) 



The interaction force in Eq. (|157|) is identical to Eq. (16) 
from Ref . , except that we have here added the terms 
quadratic in T. This equation is thus valid for any ma- 
terial as long as {d, Xt p } 3> R- In the limit where the 
separation is the largest scale, i.e. for d {R, \t p }, the 
evanescent contribution in Eq. (|157[) decays like d~ 3 [27| 
(compare also Ref. [16j for the situation of an atom and 
a plate). The self force becomes 
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(160) 

This is identical to Eq. (19) from Ref. [27|, again up to the 
terms quadratic in T ■ F$ behaves in this limit similarly 
as F 2 in Eq. (|148|l . i.e., it can oscillate for dielectrics as 
a function of d, falling off at large separations as 1 jd. 

We now repeat the expansion of the forces for small 
temperatures and dielectric spheres given in Ref. [13] ■ 
For a dielectric sphere and plate, we can employ 
Eqs. (|149[) and (|150[) . to obtain the leading behavior at 
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low temperatures (At p 3> {Xo,R}, but not necessarily 
At p ~^>d). The d independent part now becomes, 
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Fp ev can be analyzed in the following two limits, corre- 
sponding to expansions of the function f ev , 
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In the opposite limit, with At p 3> {d, Rj, Ao}, we have 
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(163) 



Equation (|T!)2"j) is similar to Eq. (12) in Ref. [l6|. As it 

is the case for F 2 , to leading order in temperature the 
self part Fg does not oscillate. For d ^> At s 3> {-R, Ao}, 
we have F s ^ cx 1/d 6 , while for X Ts > {d, A } 
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which is identical to Eq. (|163|) when interchanging real 
and imaginary parts for r p and a. 



XII. APPLICATIONS: STRONG 
NON-EQUILIBRIUM FORCES FOR 
NANOSPHERES 



Sections IXIBI and IX Bl give asymptotic formulae for 
forces and transfer for the sphere-plate geometry. Here, 
we present two examples, focussing on cases where the 
total forces allow for stable levitation when including the 
gravitational force acting on the sphere. In all cases, 
{d, At} S> R is assumed, such that we can use the corre- 
sponding approximations. The equilibrium force at finite 
temperature is computed from Eq. (17) in Ref. ffl\ . 



A. Metal sphere in a hot environment 

We start with a metal sphere, choosing aluminum for 
its low density, in front of a SiC plate. The Drudc 
model of Eq. (|128[) is used to model aluminum, with 
uj p = 12.04 eV and u T = 12.87 x 10" 2 eV [z|. For 
the dielectric response of SiC we use, [79j 
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where = 6.7, u LO = 0.12 eV, u TO = 0.098 eV, 
7 = 5.88 x 10 -4 eV. Both functions have a sufficient range 
of validity for our purposes. We find that the distance 



independent part of the interaction force in Eq. (| 157[) be- 
comes comparable to the weight of the sphere (aluminum 
has a mass density of 2.7g/cm 3 ), at roughly T = 2700 K 
(for R = 90 nm) . In Fig. [5j we show that when gravity 
is included a stable point of zero force, i.e., a levitation 
point, is possible in this regime. The main part of the 
figure shows the net force on the sphere hanging below 
the plate, at T p = 300 K and T cnv = 2862 K. In this 

situation, Fp S ^ enters Eq. (|77p with a minus sign, such 
that it becomes attractive for d — > oo, and thus almost 
balances the gravitational force (which in this setup cor- 
responds to a repulsive force). The near field part of Fp^ 
[see Eqs. (|159[) or ()163|) ] is now positive, leading to the 
repulsive barrier seen in the figure. At the shortest sep- 
arations d, the attractive equilibrium part is dominant. 

The gravitational force can be fine-tuned by using a 
spherical shell, which has (almost) identical optical prop- 
erties to a solid sphere as long as the shell thickness is 
large compared to the skin depth of aluminum (around 
20 nm) . Figure [5] shows the situation for outer and in- 
ner radii of 73 and 23 nm, respectively. To illustrate a 
potential experiment, we also show the periodic motion 
of the sphere (the solution to Newton's equation of mo- 
tion), starting from d(t = 0) = —4 fim with zero initial 
velocity in the inset of Fig. [5] The sphere oscillates in the 
potential minimum on a timescale of tens of milliseconds. 
Finally, we note that the temperature of the sphere plays 
almost no role, as the self force in Eq. f|160[) is negligible 
due to the small emissivity of a metal sphere (see Fig.EJ. 
This is of advantage, as the temperature of the sphere 
cannot be controlled in such a situation. We note how- 
ever that due to the proximity to the plate, it is likely 
to attain a temperature close to the plate temperature of 
300 K, which avoids potential melting of the sphere. 



B. Bouncing hot dielectric sphere over a room 
temperature plate 

In Ref. [29[ it was shown that dielectric materials 
can support strong non-equilibrium near field forces 
when their resonances are slightly detuned (see also 
Refs. [H, HO]). Here we demonstrate that this can lead 
to a repulsive non-equilibrium force on a small sphere 
that can exceed gravitational forces at moderate temper- 
atures. We use for sphere and plate the following oscil- 
lator models, with one resonance each in both infrared 
and the optical range, 
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The parameters used are given in Table HI and resem- 
ble realistic values [8(|. Most importantly, we have de- 
tuned the infrared resonances, which strongly changes 
the non-equilibrium forces (we found the strongest ef- 
fects for a detuning by 1.19, see Table [TT]) . The inset of 
Fig. [5] shows the total (Casimir and gravity) force for 
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FIG. 5: The net force (including gravity ) on a spherical 
aluminum shell (outer and inner radii of 73 and 23 nm) 
below a SiC plate, for T p = 300 K and T cnv = 2862 K 
(•^T onv = 800 nm). (T 3 plays no visible role.) The force is 
normalized by the gravitational force Fa, with positive val- 
ues pointing up towards the plate (i.e. gravity gives a negative 
contribution). The point of zero force at smaller d is stable; 
the inset shows oscillations around this point when the sphere 
is initially put at d = —4 (im with zero velocity. 



a sphere of R = 60 nm above a plate (gravity being 
attractive) for temperatures T p = T cnv = 300 K and 
T s = 916 K. The mass density of the sphere is assumed 
to be 2g/cm 3 . The total force is equal to the gravita- 
tional force at d — > oo (because all Casimir contributions 
vanish asymptotically in this case), and becomes repul- 
sive for separations below one micron due to the non- 
equilibrium force Fg S ^ (which is stronger than the gravi- 
tational force). For comparison, the dashed line gives the 
total force for T p = T cnv = T s = 300 K, which is purely 
attractive. We emphasize that the non-equilibrium re- 
pulsion is not due to radiation pressure but arises from 
near field (evanescent) effects. 

Again, mimicking a potential experimental realization, 
we show in the main part of Fig. [6] the solution to New- 
ton's equation of motion for a sphere dropped from a 
height of d(t = 0) = 800 nm. The falling sphere bounces 
back due the repulsive barrier, oscillating on a timescale 
of milliseconds. For comparison, we include the trajec- 
tory in thermal equilibrium, where the sphere just drops 
onto the plate. Of course the sphere will lose energy due 
to heat radiation and transfer and cool off in time. In- 
terestingly, the (near field) transfer to the plate [where 
we can use Eq. (|136l) with i max =l] is also a strong func- 
tion of the detuning parameter, and is maximal at around 
the value of 1.19 (making the transfer to the environment 
negligible in this case) . We have solved the coupled equa- 
tions for the time dependent trajectory and temperature. 
Suppressing dependencies on time independent parame- 
ters T p and T env , the two equations are 

Ts = _?^mm t (167) 
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TABLE II: Parameters for the oscillator model of the dielec- 
tric functions of the sphere and plate. uj a , Q a , -y a and T a are 
given in rad/sec, the remaining parameters are dimensionless. 
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FIG. 6: Inset: The net force on a dielectric sphere of radius 
60 nm above a dielectric plate, including the gravitational 
force Fg on the sphere (which is negative), for T v = T env = 
300 K and T 3 — 916 K. The strong non-equilibrium repulsion 
is due to near field effects, not radiation pressure, and leads 
to a stable zero force point (including gravity). The dashed 
line shows the total force for T p = T cnv = T s = 300 K. The 
main figure shows the trajectory of the sphere, when dropped 
from 800 nm. Points show the trajectory including cooling 
down of the sphere, where the length of the vertical bars is 
proportional to T 3 (t) — 300 K in arbitrary units. 

mJ(t)=F{T s {t);d(t)), (168) 

where m s is the mass of the sphere. The heat capacity 
of the sphere is estimated from n s /m s = 800 J/(kg K), 
which is a realistic value for solid materials. The solution 
in Fig. [6] indicates that the sphere falls slightly further 
due to cooling down, but still bounces back before even- 
tually cooling too much to counterbalance gravity. The 
vertical bars denote the difference T s — 300K in arbi- 
trary units. The distance dependence of the cooling rate 
(which also depends on T s ) is hardly visible. We note 
that the time to eventual drop of the sphere is roughly 
four times longer for the hot sphere (T s (t = 0) = 916 K) 
compared to T s (t = 0) = 300 K, making such an experi- 
ment sensitive to non-equilibrium effects. 

While the examples presented in this section may prove 
demanding to actual experimentation, we believe that 
they provide valuable new insights into the physics of 
non-equilibrium fluctuations. Also, compared to the ex- 
amples in Ref. j27[, we have increased the ratio of non- 
equilibrium repulsion to gravity by almost six orders of 
magnitude, and it is well possible that it can be increased 
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even further by choosing appropriate materials or com- 
pounds. 
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Appendix A: Equilibrium force formula from the 
field correlator 

The trace formulae for the non-equilibrium force in 
Eqs. ([74")) and (|75"|) were derived starting from Eq. (fT3"|). 
This starting point can be further supported by show- 
ing that it leads to the known result for the equilibrium 
Casimir force, as demonstrated in this Appendix for two 
objects (1 and 2). In equilibrium, the field correlator is 
given by Eq. Q , and all we need is the Green's function 
of the system. It is found by starting from object 1 in 
isolation with Gi = (1 + GoTi)Go (see Eq. (|2"Tj) ). and 
inserting object 2 by use of the operator Oi in Eq. (j4"T|). 
as 



G = (1 + G T 2 ) 



1 



1 - G0T1G0T: 



■(l + G Ti)G . (Al) 



In applying Eq. (|73|) . we note that we may as well con- 
sider the complex conjugate of the integrand, in which 
case we have G^" 1 , a better fit to Eq. (|A1[) . Furthermore, 
we note that by symmetry, only terms with an even num- 
ber of T matrices contribute to the force. Using the ar- 
guments given below Eq. (f73"T) . we find the forces on the 
two objects to be (note that the field correlator in Eq. ((4]) 
is real), 



p(l,eq) 



■in 



did 



-p(2,eq) 



Im Tr <! VG T 2 



2h f 00 , 

= — / duj 

I" JO 



1 

T + 5 
1 



Noting that V r Go(r,r') = — V r <Go(r, r'), one can see 
that subtracting Eq. (|A3[) from Eq. (|A2[) yields a result 
which can be written as a derivative with respect to the 
position of object 1, as 



P (i,e,) j_ r duJ 



1 



1 



ImTrlV^logll-GoTiGoTa]}. (A5) 
This derivative can be taken out of the trace, leading to 



V 0l T =V , - 



or 



1 



1 



ImTr {log [1 - G TiG T 2 ]} . (A6) 



Here we have introduced the Casimir free energy T Q, 
and the result in Eq. (|A6|) is consistent with previous 
work. Using Eqs. (f99|) and (1821) . wc find the even more 
familiar result in terms of partial waves, 



7T 



die 



(A7) 



x ImTrlog[Z-W 21 TiW 12 T2]. 

Appendix B: Plane waves basis 
1. Free Green's function 



The plane waves basis is most conveniently chosen with 
the symmetry axis pointing along the z-direction, for 
plates lying in the xy plane; xi and kx are the spa- 
tial coordinate and the wave- vector perpendicular to this 
axis, respectively. The plane wave basis is fundamentally 
different from coordinates with a singular point (such as 
the spherical basis) and we give two representations of 
the Green's function with different properties. Wc con- 
sider the following vector eigenfunctions, 



1 - G TiG T 2 
1 



M±Jx ± ,z) = 
GoTi \ , (A2) N±Jx x ,z) 



2^|k 



{kk y - yk x ) e l 



2Vfe|k 



(zLxk x k z i ykyk z + z/cj 



z -L 



1 



(Bl) 



(B2) 



ImTr VG TiG T 2 



1 



1 - G TiG T 2 



(A3) 



where k = (kj_, k z ) T , with k z = <J ^ — k\. We consider 

two sets of eigenfunctions, for reasons discussed below, 
and start with elementary right and left moving waves, 



In Eq. (|A3|) . we have re-summed the expansion of the in- 
verse operator. The above two equations satisfy F( 1,e? ) = 
„jr , (2,e9) ag expected, which can be seen by expanding the 
inverse operators. We can hence write 



Ek e ,Wr) = P k > ± ,z), 



p(l,e 9 ) 



e i) — jr , ( 2 ' e 9) 



(A4) 



< 



(B3) 
(B4) 

(B5) 
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z > 

2E™s (r) z<0 



(B6) 



E 



s.Pk 



V2 



Wp.kJrJ + sE^Jr)]. (B8) 



We also define waves of definite parity under reflections 
at the z = plane, that carry an index s = ±, 



l 2 



.Pk I 



(r) = — p- E^p k (r) + sEjpjj (r) , (B7) The free Green's function can now be written as 



G (r,r') 
Go(r,r') 



' E E 

P=N,M j=L,R 

' E E 

P=N,M s=± 



d 2 k ± j E^, k Jr) ® Ef« _ k Jr') |z| > |z'|, 
(2tt) 2 \ E^_ k Jr)®E°£ k Jr') |z| < |z'|, 

rf 2 fc± f E°£ k Jr) E^_ k Jr') |z| > |z'|, 
(2tt)2 \ E r ^_ k Jr)®E°£ k Jr') \z\ < \z>\. 



(B9) 
(BIO) 



where L — R and P = L in the first line. We note that 
both representations are of the form of Eq. ([75)1 , and obey 
Eq. The first variant is convenient for study the 

problems involving only one side of a (thick) slab, as done 
in this manuscript. These waves have fi = (j, P, kjj and 
<r(fi) = (j,P, — kj_), but they do not obey the required 
properties under complex conjugation in Eqs. (|85 p — (|87 [) . 
On the other hand, the waves with definite parity ful- 
fill all properties required, i.e., Eqs. ([79j). ([gU]>. ([85]> - 
(|87|). for our final formulae for transfer and forces to 
be directly applicable. They have /i = (s,P, kj_), and 
cr(/i) = (s, P, — kj_), and their phases are 

= [0(w/c - |kj_|) - * (-l)Ci-«)/2 (|k x | - w/c)]. 

(Bll) 



2. Fresnel coefficients 

The Fresnel reflection coefficients r p for reflection are 
given in Ref . |(36| , 



r M (fcx,w) = 



/' 



0, and <A, are 



Mlm 



Eleg 



ii -r V x ry, 



(CI) 



VxrK 



E£f = -V x E™e , 

JVtm , , Mini 

EOUt C T-7 TTlOUt 

JVZm — —V X & Mim- 



(C2) 
(C3) 

(C4) 



ji is the spherical Bcsscl function of order Z, and /i; is 
the spherical Hankcl function of the first kind of order I. 
Y l m (9,(f>) are the spherical harmonics, where we use the 
standard definition according to Ref. [(36| (which is dif- 
ferent from the one used in Ref. [7(|). The free Green's 
function Go is then given by Eq. ([79)1 with /i = {P, I, m} 
and a(p) = {P,l,-m}, and £ M -> Ep Ez=i EL=-r 
These definitions fulfill Eqs. (JHfi) and flSSHETJ)- The 

(B12) 

phase function is unity throughout, e l< ^' = 1. 



r N is obtained from r M by interchanging /i and e. In 
Sec. IIX Al we give the relation between the Fresnel coef- 
ficients and the T matrix elements of a plate. 



2. T-matrix of a sphere 



Appendix C: Spherical Basis 

1. Partial waves and free Green's function 

Here we adopt a wave expansion similar to Ref. [70j . 
where the waves, depending on spherical coordinates r, 



The matrix elements of T for a sphere of radius R, 
as defined in Eq. (|82[). are well known [7{J and some- 
times referred to as Mic coefficients. Considering for 
simplicity spheres with isotropic and local e and (i, ren- 
ders the matrix T v P l ^ m diagonal and independent of m, 



n-P'P 
' Vim' ra 



7~i F Spp'8u'5 mm ' . The matrix element can 
be conveniently written in terms of R* = Ruj/c and 
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/sJIRuj/c, as 





[R*ji(R*)] 




~R*ji(R*) 




[K*hi(R*)] 




~R*ji(R*) 








(C5) 



77 M 



77 follows from 77/ by interchanging and e. 

Appendix D: Cylindrical Basis 
1. Free Green's function 

The cylindrical coordinates are denoted by p (radial), 
z (along the cylinder axis), and azimuthal angle <$. The 
free Green's function is now expanded in terms of cylin- 
drical waves, which, using a notation similar to the one 
in Ref. [70[, are defined by 



E 



M. 



-prog 

Eout 



-V X KJ, 



U) 
c , 



= - V x E out 



M,n,k I ' 
out 

M.n,k z ) 



(Dl) 

(D2) 
(D3) 

(D4) 



with k„ 



- — kg, and 

=Uk eQ y k > z+m t (D5) 

(q, z, 0) =H « (fc^e**-^*** (D6) 

J„ is the Besscl function of order n, and £Zn is the 
Hankel function of the first kind of order n. The free 
Green's function Go is then given by Eq. ([79]) with 
fx = {P,n, k z } and = {P,-n,-k z }, and ~> 

EpE^=-oo/lF- These dcf initions fulfill Eqs. (Sty 
and J5S ) -l|g7 ) with the phase 

e**" = [6(w/c- |*,|) +(-l)( n+1 )9(|fc z |- W / c )], (D7) 

where the 8-function is used to stress that the phase is 
unity for propagating waves. 



2. T-Matrix of a cylinder 

The T-matrix for an infinite homogeneous cylinder is 
diagonal in n and k z , but not in polarization P [52|. 
The coefficients T^ k p take a lengthy form and are not 
reproduced here. They can be found in Refs. [52l. l55l. l75j : 
in Ref. (55j the matrix elements are given in precisely the 
notation used here. Ref. [55J also gives T p k p for uniaxial 
materials. 



Appendix E: Conversion matrices 

1. Plane waves to spherical waves 

The outgoing plane wave eigenfunctions in Eq. (|B5 
are expanded in spherical waves in the following way, 



7 E° fl V(x 1 ,z)= J2 D lmP'P^ph m , (El) 

P',l,m 



with the matrix elements expressed in terms of Legendre 
Polynomials P z m , as 
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lmMMk.j 



' 4vr(2Z + 1)(Z - m)\ IT\k ± \ 



(E2a) 



ImNMk 



yPlFV Z(Z + l)(/ + m)! c|k x | 



AmMJVk , = Di mN M\t I • 



(E2b) 

(E2c) 
(E2d) 



is the the angle of with respect to the x-axis. 

2. Spherical waves to spherical waves 

Outgoing spherical waves can be expanded in regular 
spherical waves with respect to a different origin, shifted 
by d. In the cases considered in this paper, the transla- 
tion can always be chosen along the z-axis, such that 

Ep? m (r) = £W± Pi ^d)E r ;s, m (r ± dz). (E3) 

P'V 

For example, if the coordinate system of object 2 is cen- 
tered at £>2 = Oi - dz as considered in Eqs. (|14ip . (|142p . 
(fT45]) and JUS), one has W 21 = W+ in Eq. J99). The 
elements of W are [7l| , 



U ± 

"P'PA'lm 



E 



i(f + l) + i'(f' + l)- 1/(1/ + 1) 



^fimd— (1 — dppi) 
c 



(E4) 



with the function 

4,»(^("i) m ^ !,± T^i)i 



f (2f + 1)(2Z / + 1) 
J(J + 1)J'(Z' + 1) 
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Eq. (una), 



- Pz-.P'Pi'im = — - Sp>p) &i>,i a(l, m) 
c 

+ 5p' P [-b(l, m) Si>j + i + b(l', m) 5i>+i,i]} 
useful relation for the matrix — p z = ddV(di)\d=o (see with a(l,m) and b(l,m) defined in Eqs. (|143[) and (|144l) 



The regular part V, as used in Eq. (fHTjl . is obtained from + Sp>p[-b(l, m) 8 v ,i+i + b(l', m) Sy+ij}} , (E6) 

U + by replacing h v in Eq. (|E5|) with j„ [rij . We find the 
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We note that Eq. (|135l) has an additional factor of 8n 
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